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ABSTRACT 

We present a multi-wavelength study of three star forming regions, spanning the age range 
1-14 Myr, located between the 30 Doradus complex and supernova SN1987A in the Large 
Magellanic Cloud (LMC). We reliably identify about 1000 pre-main sequence (PMS) star 
candidates actively undergoing mass accretion and estimate their stellar properties and mass 
accretion rate (M). Our measurements represent the largest M dataset of low-metallicity stars 
presented so far. As such, they offer a unique opportunity to study on a statistical basis the 
mass accretion process in the LMC and, more in general, the evolution of the mass accretion 
process around low-metallicity stars. We find that the typical M of PMS stars in the LMC 
is higher than for galactic PMS stars of the same mass, independently of their age. Taking 
into account the caveats of isochronal age and M estimates, the difference in M between 
the LMC and our Galaxy appears to be about an order of magnitude. We review the main 
mechanisms of disk dispersal and find indications that typically higher M are to be expected 
in low-metallicity environments. However, many issues of this scenario need to be clarified 
by future observations and modeling. We also find that, in the mass range 1-2 M Q , the M of 
PMS stars in the LMC increases with stellar mass as M ocM* with foal, i.e. slower than the 
second power low found for galactic PMS stars in the same mass regime. 

Key words: Galaxies: Magellanic Clouds - stars: formation - stars: pre-main-sequence - stars: 
circumstellar matter 



1 INTRODUCTION 

The mass accretion rate (M) of young stars is a key parameter to 
constrain the models of both star and planet formation. M is an 
indicator of the presence of gas, and hence dust, in the inner cir- 
cumstellar disk and, as such, it is affected by the disk structure and 
evol ution as well as by the formation and migration of planets (see, 
e.g. JCaivet et al. 2000). Thus, it is of particular interest to deter- 
mine the evolution of M as a star approaches its main sequence, 
whether M and stellar mass are correlated, and whether it is af- 
fected by the specific conditions of a given star formation event, 
such as the chemical composition and density of the parent molec- 
ular cloud, by the proximity of hot massive stars, etc. The measure 
of M generally relies on the study of continuum veiling, UV excess 
emission or profile/intensity of emission lines (usually Ho-, Pa/? or 
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Bry), which requires medium to high resolution spectroscopy of 
individual young objects. As a result, the data collected so far are 
limited to nearby (d<l-2kpc) star forming regions. These data indi- 
cate that M is a decreasing function of s tellar age (Muzerolle et al. 
l200d : ISicilia-Aguilar et alj|2005l l2006t |20ld : iFedele et alj|20ldx 
in lin e with the expected evolution of viscous discs dHartmann et al.l 
1998), and is roughly proportional to the second power of the stellar 
mass dMuzerolle et alj2005l : lNatta et al]200fj| : ISicilia-Aguilar et all 
2006). However, the detection limits and the difficulty of measur- 
ing very small M do not allow us to rule out d etection/selection 
thresh olds as responsible of this trend (see, e.g-. IClarke & P ringle 
2006). Moreover, the spread of the M data exceeds 2 dex at any 
given age. Several facets of the mass accretion process appear to 
contribute to this spread: i) at any age, disks in different evolu- 
tionary stages coexist in the same star forming region (see, e.g., 
lHarvev et alj2007l : lAlcala et alj2008l : lMerin et al.l2008h : ii) M de- 
pends on the disk mass and, hence, on the mass of the forming sta r 
(Whit e & Hillenb rand 2004; lAlcala et alj|2008l : [kirn e t al. 2009); 
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iii) the accretion process in young stars, and the accretion-related 
emission lines, are subject to variations on a timescale of a few 
hours to years or even decades dHerbstll 19861: lHartigan et al"]|l99ll : 
iNguven et al]|2009h . Another limitation of the galactic objects for 
which M measurements are curr ently available is that they have 
essentially solar metallicity (e.g. Padgetti ll996h , because they are 
located in nearby star forming regions. 

In the past two decades, many authors have demonstrated the 
potential of deep Hubble Space Telescop e (HST) optical imag- 
ing for star formation studies. Recently. |Pe Marchi et al. 
(hereafter Paper I) presented a novel method to identify pre-main 
sequence (PMS) objects actively undergoing mass accretion and 
determine their M. This method relies solely on HST imaging, 
which can be easily obtained for a large number of PMS stars at 
once and for distant star forming regions, well beyond a few kpc 
and with different metallicity, allowing us to overcome some of 
the limitations affecting techniques based on spectroscopy. This 
method has been already adapted to the Wide Field Camera (WFC) 
of the Isaac Newton Telescope (INT) and is being used within 
the frame of INT Photometric H-Alpha Survey (IPHAS), a 1800 
deg 2 survey of the Northern Galactic Plane aiming at reliably se- 
lect classical T Tauri s tars and constrain their mass accretion rates 
jBarentsen et alj201 ll) . We have now started a pilot project aiming 
at applying this technique to the existing and extensive HST pho- 
tometry of star-forming regions in the Milky Way (MW) and the 
MCs. In this paper we present the results of this study in the Large 
Magellanic Cloud (LMC). 

Because of it r elative proximity (51.4+1.2 kpc; Panagia 1991 ; 
Panagi a etaljl99lh and low-metal content (Z=0.007; Mae der et all 
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19991 and references therein), the LMC is a very interesting test- 



case for star formation studies and, specifically, for mass ac- 
cretion studies, because it offers a unique opportunity to study 
a resolved accreting stellar population in an environment with 
a chemical composition different from the MW (Z«0.018; see 
lEsteban & Peimbertl 19951) . 

Investigating the M vs. metallicity dependency is of special 
interest for several reasons: 

(i) The time-scale of the mass accretion process is strongly con- 
nected with the disk lifetime, which h as been poorly stu died in 
low-metallicity environments (see, e.g.. lYasui et al.l 12009). Thus, 
studying the M evolution might help clarifying how different disk 
dispersal mechanisms act at lower-metallicity. 

(ii) Low-metallicity disks can be used as a benchmark for 
evolved solar-metallicity dis ks, where dust grains ma y have suf- 
fered strong coagulation (see lErcolano & Clarkel201cL and discus- 
sion at the Ringberg Accretion Workshop 201 10). 

(iii) The evolution of M is affected by possible planet formation 
in the disk and, hence, it might provide new important clues on 
on the so called "planet-metallicity correlation". The probability 
of a s tar hosting a planet appears to increase with stellar metal - 
licity dGonzalej|l997l ; lFischer & ValentfcOollSantos et al]2003l) . 
However, a few planets have been found around metal -poor stars 
(see, e.g. ISantos etai1l2007l : ISetiawan et al]|2003l. 12010) and the 
frequency of planets around metal-poor stars compared to solar- 
metallicity stars is still uncertain. Thus, the "planet-metallicity cor- 
relation" could be a simple selection effect. Indeed, surveys for 
extra-sol ar planets are mainly based on the radial ve locity (RV) 
method dMavor & Ouelozlll995l ; rMarcv & Butlerll 19981 . and refer- 
ences therein) and RV variations are much harder to detect in stel- 



lar spectra presenting few metallic lines. As a consequence, only 
a few planets have been found around very metal poor stars and, 
because of this indication, metal-poor stars are systematically ex- 
cluded from RV surveys. 

This paper is organized as follows: Sect.|2]describes the selec- 
tion of the LMC fields among those available in the HST archive 
and the data reduction. In Sect. [3] we apply the method described 
in Paper I to select PMS objects and calculate their M. The results 
are presented in Sect.[4]-[8] where in particular we discuss the mass 
dependency of M, its evolution and metallicity dependency in rela- 
tion with disk dispersal mechanisms such as planet formation and 
disk photo-evaporation. The conclusions of this work are given in 
Sect.[9] 



2 DATA SAMPLE 

We selected three 162"x 162" areas in the LMC located between the 
30 Doradus star forming complex and supernova SN1987A (Fig- 
ure[TJ. All fields were selected from the HST Archival Pure Paral- 
lel Project; the proposal ID and Principal Investigator (PI) for each 
dataset are reported in Table [T] The observations were performed 
from February 1999 to November 2003 with the Wide Field Plane- 
tary Camera 2 (WFPC2) of the HST using the F606W (V), F656N 
(Ho-) and F814W (I) filters. The three fields were selected following 
two basic criteria: i) exposure time long enough to allow the detec- 
tion of objects down to 0.8-1 M G at the distance of LMC in the 
three filters; ii) different environmental conditions, such as vicinity 
to OB stars and amount of interstellar gas and dust, in order to study 
the behavior of M in star-forming regions with different properties. 

The basic data set used for this paper was retrieved from the 
HST archive at the Canadian Astronomy Data Centre (CADC). 
The sources identification and photometry was performed on the 
stacke d images using the SExtractor 2.5 tool bv lBertin & Arnoutsl 
dl996h . SExtractor exploits the aperture photometry technique, 
which is the fastest and best approach for not-crowded fields as 
those considered here, and the background is locally estimated in 
an annulus surrounding the star. Following the prescriptions o f 
the HST Data Handbook for the WFPC2 dBaggett et all l2002h . 
we adopted a 0.5 arcsec aperture radius for magnitude extraction 
and then calibrated the instrumental magnitudes to the VEGA- 
MAG photometric system using the zeropoints listed in Table 5.1 
of the handbook. The SExtractor morphological parameter "extrac- 
tion flag" (FLAGS) was used to clean the catalogs from spurious 
detections, such as saturated or truncated sources too close to the 
image boundaries, etc. Table [T] summarizes, for each of the three 
fields, the total exposure time, the saturation limit and the limiting 
magnitude achieved at the lOcr level in each filter. The exposure 
times are not uniform across our dataset, producing different satu- 
ration limits and limiting magnitudes for each filter/region. Consid- 
ering the distance to the LMC (~50 kpc) and the typical age of our 
re gions (Tabled and u sing the evolutionary models for PMS stars 
bvlTognelli etal.U201lh for the metallicity of the LMC (Z=0.007; 
Maede r etalj|l999l" and references therein), we concluded that our 
observations are complete above 1 M Q for the three regions in all 
the filters. Thus, in Sect. |5l6| we will limit our analysis to the mass 
range 1-2 M G . 

In addition to the above fields, we have used WFPC2 
photometry of a field around supernova SN1987A studied by 
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See http://archive.stsci.edu/prepds/appp for details. 
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Figure 1. R-band image of a 40'x40' area in the LMC form the Sloan Digital Sky Survey. The polygons display the four regions observed with HST/WFPC2 
analyzed in this paper. The position of 30 Doradus star-burst cluster and the supernova SN1987A are also indicated. 



Romaniello dl998h . IPanagia et all l l2000h and iRomaniello et al.l 
1 20021) . These observations are centered at R.A.=5 ft 35 m 28 s and 
Dec.=-69' ,< ' s 16'" 12 s and cover a field of 9.16 square arcmin. The 
mass accretion rate for the PMS stars in this field were already de- 
rived in Paper I. 



3 SELECTION OF PMS STARS 

We select PMS stars in each field and determine their mass accre- 
tion rates using narrow-band Ho- (f656n) and broad-band V (f606w) 
and I (f814w) dereddend photometry. In the following sections we 
describe the main steps of our selection procedure. For a more de- 
tailed description of this method we refer the reader to Paper I. 



3.1 Reddening 

From a visual inspection of the images, we find evidence of patchy 
nebulosity in our fields and we therefore can expect a certain 
amount of differential reddening. 

To test this hypothesis, we divided each field in 22 regions, 
each 32"x32" in size, and obtain the V vs. V-I colour-magnitude 
diagram (CMD) for each sub-image. We used the zero age main 
sequence (ZAMS) bv lMarigo et al.l J2008D for the metallicity of the 
LMC as the reference unreddened CMD (Figure |2j. The ZAMS 
was shifted to the distance of the LMC and reddened by the amount 
corresponding to the interven ing absorption along the line of sight 
due to the Milky Way, which Fitzpatrick & Savage! d 1984b quanti- 
fied in A v =0.22, in turn corresponding to E(V ? I) =0.1. 

Within each sub-image we derived for all stars a single value 
the reddening relative to the ZAMS in the following way. For each 
star with 18<V<24 and 0<V-I<1.3, we calculated the distance from 
the ZAMS along the reddening vector, defined by the extinction 
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Field 


HST proposal ID 


PI 


RAJ2000« B(er 
(hh:mm:ss) 


DECJ2000„„,„ 
(dd:mm:ss) 


Filter 


(min) 


Mag Sat. 


Mag lOcr 


#1 


8059 


Stefano Casertano 


05:34:34.03 


-69:14:01.50 


F606W 
F656N 
F814W 


17 

73 
39 


17.00 
15.00 
17.00 


24.30 
21.00 
23.80 


#2 


9634 


James Rhoads 


05:35:34.93 


-69:07:32.36 


F606W 
F656N 
F814W 


49 
68 
35 


17.50 
15.50 
16.00 


25.20 
21.30 
24.30 


#3 


8059 


Stefano Casertano 


05:39:40.57 


-69:10:48.43 


F606W 
F656N 
F814W 


8 

57 
15 


16.00 
14.30 
15.50 


23.80 
20.00 
23.50 



Table 1. HST proposal ID, PI name, total exposure time, saturation limit and limiting magnitude at 10cr level for the three HST-WFPC2 pointings used in this 
paper. 



law as measured specifically in the LMC bv lScuderiet~ai1 i ll996h . 
Each one of these distances is the resultant of two components, 
namely E(V-I) along the abscissa and Ay along the ordinates. We 
then looked at the distribution of Ay values inside each cell and 
took the 17th percentile as a conservative lower limit to the value 
of Ay for the whole cell (see example in Figure |3j- If the redden- 
ing distribution were Gaussian, this would correspond to lcr. The 
purpose of this approach is to minimize the number of stars that 
after the correction become bluer than the ZAMS and to keep them 
to a fraction compatible with our photometric uncertainties. This 
choice ensures that we are not overestimating the reddening cor- 
rection for the majority of the stars in the field. It also avoids an 
improper overestimation of the ages due to an excessive reddening 
correction (see reddening lines in Figur^2}. 

The relative Ay values of the sub-images span a range of 
0.4 mag in our field #1, 0.6 mag in our field #2 and 1.2 mag in our 
field #3. Thus, differential reddening is more significant in field #3, 
which is indeed the closest one to the 30 Doradus star-burst cluster 
and is still associated with conspicuous gaseous/dusty nebulosity 
(Figure |9j- 

We correct the magnitudes of each source in each filter assum- 
ing the Ay estimated for the relative sub-image and the extinction 
law bv lScuderiet"aT1 < ll996h . This statistical method, though pro- 
viding only a rough estimate of the visual extinction toward the 
given source, is suitable for our purposes. Indeed, in Paper I we 
have demonstrated that for an average LMC star-forming region 
omitting the extinction correction would result in a 10% overes- 
timate of the mass accretion rate. This error, albeit systematic, is 
smaller than most other systematic uncertainties and comparable to 
the measurement errors. 



3.2 Selection of stars with Ha excess emission 

We select PMS star s on the basis of their Ha emission line 
(Whit e & Basrl| [2003). The traditional approach to search for Ha 
emitters is based on the use of the R-band magnitude as an indica- 
tor of the level of the photospheric continuum near the Ho- line, so 
that stars with strong Ha emission will have a large R - Ha color. 
Since we do not dispose of R-band photometry for our fields, we 
use measurements in the neighboring V and I bands, following the 
approach described in Paper I. 

Figure|4](top panel) shows, as an example, the V-Ha vs. V-I 
color-color diagram for our field #1. We use the median V - Ha 




12 3 4 
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Figure 3. Example of Ay distribution for stars in a sub-image of our field 
#1. The dashed line indicates the 17th percentile of the distribution. 

dereddened color of stars with small photometric errors in the V, I 
and Ha bands (S 3 = ^(&V 2 + AI 2 + AHa 2 )/3 < 0.15) as a func- 
tion of V - I to define the reference template with respect to which 
the excess Ha emission is sought (line in Figure|4j. Since the con- 
tribution of the Ha line to the V magnitude is completely negligi- 
ble, the magnitude AHa corresponding to the excess emission is 
simply: 

AHa = (V- Ha) obs - ( V - Ha) Icf ( 1 ) 

where the superscript "obs" refers to the observations and 
"ref" to the reference template. We select a first sample of stars 
with excess Ha emission by considering all those with <5 3 <0.15 
and V-Ha color above the reference template by 5 times their in- 
trinsic photometric uncertainty (AHa > 5 x VAV 2 + AHa 2 ). Note 
that this threshold is somehow arbitrary; it affects the total number 
of selected PMS star candidates but not the typical M of the se- 
lected PMS population, on which the conclusion of this paper are 
based (Sect.[6]>. 

Then, we derive the equivalent width of the Ha emission line 
(EW//„) from the measured AHa as: 

EW Ha = RW x [1 - KT 04xA/,ff ] (2) 

where RW is the rectangular width of the filter, which depends on 
the characteristics of the filter (see Table 4 in Paper I). We also 
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Figure 2. Example V vs. V - I CMDs for one of our WFPC2 fields (field #1) divided in to 22 32"x32" regions to study the spatial variation of the interstellar 
extinction over the observed field. The solid line is the ZAMS bv lMarigoetal] {2008). The dashed lines show the direction of the reddening vector and, 
together with the dashed dash-dotted lines, define the region of the CMD uncontaminated by stars in the red clump. The visual extinction assigned to each 
region is reported in the top-left of each panel. 



derive the effective temperature (T e ff), radius (R*) and luminosity 
(L*) from the dereddened V - I color and V magnitude. We adopt 
the di stance of 51.4 kpc to the LMC jPanagia|[l99ll : IPanagia et all 
and use the tabulation of HST/WFPC2 colors as a function 
of T e// bv lFessel et al.1 < fl998l) . 

Following the prescription of Whi te & Basril d2003h and 
adopting the effective tem perature vs. spectral type scale by 
iKenvon & Hartmannl < fl995h . we finally consider as probable PMS 
stars those objects with (Figure [4] lower panel) EW# ff >3 A for 
stars earlier than K5 (i.e., T e// >4350 K), EW Ha >10 A for K5-M3 
stars (340CK T eff <4350 K) and EW Ha >20 A for M3-M6 stars 
(3000< T e ff <3400 K). Our photometric limits are such that we 
do not detect LMC members later than M6 in any of the fields (see 
Sect.[2jl. We also reject stars with V-I <0. These two criteria allow 
us to clean our sample from possible contaminants, i.e. older stars 
with chromospheric activity and Ae/Be stars, respectively. 

Note that dKe and dMe stars with enhanced coronal 



emission have EW tfa well below 20 A eve n during flares 
dWorden & Petersonl 19761: IBopp & Schmitzl 19781) . Thus, the crite- 
rion bv lWhite & B asri ( 2003) allows us to safely remove dMe stars, 
while some the dKe might still be included in the PMS sample. We 
could use a more conservative threshold of EW H „ =20 A for both 
K and M type stars, as was done in Paper II. However, while avoid- 
ing contamination from magnetically active field stars, this thresh- 
old might cause the rejec tion of some of the w eaker K0-K7 type ac- 
cretors (M < 10~ 9 M m /vn lWhite & Basril2003l) and, hence, an over- 
estimation of the typical accretion rates of our PMS populations, 
in particular at ages older than ~10 Myr. T hus, for K type stars w e 
choose to adopt the more relaxed criterion o f lWhite& Basril j2003l) . 

We identified 490 objects in our field #1, 325 in our field #2 
and 242 in our field #3 with Ha emission typical of accreting PMS 
stars. These numbers are of the same order of magnitude as the 
number of low-mass stars expected in each region on the basis of 
the number of OB stars (Sect. 18.4.21 and a typical IMF (a=2.35; 
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Figure 4. Upper panel: Example of V - Ha vs. V - I diagram of stars in 
one of our fields (field #1 ). The line is the reference template with respect to 
which the excess Ha emission is sought. Objects with a V - Ha excess are 
indicated as big dots with error bars. Lower panel: Ho- equi valent width as 
a func tion of V- / for stars in field #1. Objects satisfying the White ^k Basril 
120031) EWh„ criterion (dashed line) and redder than V - 1=0 (dot-dashed 
line) are considered PMS star candidates (big dots). The remaining objects 
(open circles) are older stars with chromospheric activity or Ae/Be stars. 



ISalpetd l 19551) . meaning that we are detecting a large fraction of 
accretors in those regions. This point is very important for the con- 
clusions of this paper. As we will see in Sect. [6] the typical mass 
accretion rates in the LMC are higher than in the MW; because our 
samples are representative of the entire PMS population in the stud- 
ied regions, we can be confident that this is not a selection effect. 

3.3 Determination of mass and age 

The mass (M t ) and age of each star are derived by comparing its 
L* and T e ff with PMS evolutionary models. Specifically, we used 



the P isa database of PMS tracks and isochrones dTognelli et al.l 
l201ll) for metallicity Z=0.007, as appropriate for th e LMC, and 
follow ed the interpolation procedure developed by iRomanieliol 
dl998h . which does not make assumptions on the properties of the 
population, such as the functional form of the IMF. On the basis of 
the measurement errors, this procedure provides the probability dis- 
tribution for each individual star to have a given value of the mass 
and age (th e method is concept ually identical to the one presented 
recently bv lDa Rio etall d2010|) ). 

Many caveats should be considered when dealing with masses 
and ages of PMS stars derived from evolutionary models. We will 
discuss them in detail in Sec. [7] because the uncertainty on these 
parameters, in particular on the isochronal age, could have a strong 
and unavoidable impact on our results. Here we stress that, al- 
though isochronal ages of individual objects are uncertain, given 
the presence of unresolved binaries and the stellar variability, ages 
of statistical samples are, in principle, reliabl e and the global ag e 
differences between regions are real (see, e.g. jMavne et al.l l 2007). 
This statistical approach is the one we use for the analysis illus- 
trated in the next sections. Moreover, there are several arguments 
supporting the reliability of our age estimates: 

(i) the median age estimated for our field #1 (13 Myr, Table®, 
the closest one to supernova SN1987A, is roughly consistent with 
the age independently e stimated for the prog enitor of SN1987A and 
nearby stars (~13 Mvn lScuderi et alj 19961) : 

(ii) Figure[5]shows the HR diagram for the PMS populations of 
the four LMC regions under analysis. There is a clear difference 
between the distribution of PMS stars belonging to the youngest 
(#3, 6-8 Myr) and the oldest (SN1987A, 13 Myr) of our regions, 
with region #3 being clearly a few Myr younger than the SN1987A 
region. The age of regions #1 and #2 is intermediate between region 
#3 and the SN1987A region and their PMS population concentrate 
between 4 and 14 Myr, where the isochrones become very close to 
each other and it is hard to derive an age difference based only on 
the visual inspection of the HR diagram. However, as reported in 
Table® the median isochronal ages of region #1 and #2 differ by a 
few Myr and with respect to region #3 and the SN1987A region as 
well, a nd the difference i s statistically significant. 

(hi) lElmegreenl ( [20r3rj|) showed that the age dispersion of a given 
stellar population depends on the size of the region. Taking into 
account the typical size of our region ( 162"x 162", i.e. about 100 pc 
at the LMC distance), we would expect an age dispersion of a few 
to 10 Myr inside each region, which is consistent with our findings 
(Table®. 



3.4 Determination of the mass accretion rate 



The mass accretion rate is derived from the free-fall equation as: 

log M(M e /yr) = -7.39 + log ^ + log ^ - log ^ (3) 



where L, acc is the energy released b y the accreti on process, directly 
proportional to the Ha luminosity (Dahm 2008, and Paper I): 

log L acc = ( 1 .72 ± 0.47) + log L(Ha) (4) 

Once the color excess AHa is determined, the Ha emission line 
luminosity, L(Ho'), can be immediately obtained from the photo- 
metric zero point and absolute sensitivity of the instrumental set-up 
and from the distance to the sources (see Sect. 3.1 in Paper I). 

The problems involved in me asuring M have been recently 
discussed by several authors (e.g.. iHerczeg & Hillenbrand! 120081 : 
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Figure 5. HR diagram for the PMS population in the four LMC regions under analysis. The dashed lines are the PMS isochrones by Tognelli et al. (201jl) for 
the metallicity of the LMC. All objects shown here have Ha excess emission and Ha equivalent width above the levels discussed in Sect. 13.21 



iBarentsen etai]|201ll ; lRigliaco et alj|201ll) and measurements us- 
ing different diagnostics might give different results by about an 
order of magnitude. The uncertainties of the Ha diagnostics and 
the specific method adopted in this work have been fully addressed 
in Paper I. The typical error on log M is of the order of 0.7-0.8 dex 
(see TableO and the main sources of this uncertainty are: 

(i) The uncertainty on the L, ICC vs. L(Ha) calibration relation 
(Equation [4]l, which is based on near-infrared spectra of a dozen 
accre ting members of the IC 348 cluster (see Figure 17 by [Damn 
2008). Due to the poor statistics, the resulting uncertainty on h acc is 
as high as a factor of 3. Moreover, the relation is mainly calibrated 
using stars with mass below 1 M Q and a handful of stars around 
2 Mo, while a gap in accreting stars is observed extending between 
0.8 and 1.8 M Q . Thus, this relation is more uncertain when applied 
to the mass regime considered in this paper (1-2 M Q ); 



(ii) Systematic uncertainties due to i) discrepancies in the 
isochrones and evolutionary tracks, ii) reddening, iii) Ha emission 
not due to the accretion process, and iv) the contribution of the neb- 
ular continuum to the colors of the stars. In Paper I we discussed 
in detail these four contributions and concluded that: i) systematic 
uncertainties of the order of 20% are to be expected for the mass 
due to differences between evolutionary models; ii) for an aver- 
age LMC star-forming region, omitting the extinction correction 
would result in a 10% overestimate of M. This error is smaller in 
our case, because our photometric catalogs are dereddened taking 
into account the effect of patchy absorption (Sect. 13. il l; iii) possible 
sources of Ha emission other than accretion are chromospheric ac- 
tivity, whose contribution is 2 orders of magnitude less than what 
we measure, Ha emission along the line of sight arising in knots 
of ionized H (e.g., very compact HII), which we removed by a vi- 
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sual inspection of all objects with excess Ha emission, and gaseous 
regions around the object ionized by an external source (e.g., hot 
ionizing stars), which, according to the calculation presented in Pa- 
per I (Table 1), is negligible in our case because there are no OB 
stars in the close vicinity (<5") of our PMS star candidates (see 
Sect. 18.4.21 ). We can, therefore, conclude that in our fields the con- 
tribution of Ho- emission generated by sources other than the ac- 
cretion is negligible; iv) if present, nebular continuum will add to 
the intrinsic continuum of the star, thereby affecting the observed 
broad-bands colors of the source. As shown in Paper I (Table 2), 
for the typical spectral range of our candidates (G-K), the effects 
of the nebular continuum on the V - I color of PMS stars remains 
insignificant even for the PMS objects with the highest EW H[ ,; 

(iii) Other possible systematic errors on the mass and age esti- 
mates due to inaccuracies in the models' input physics, which we 
will discussed in Sect. [7j 

(iv) As for the statistical uncertainty on other quantities in Equa- 
tion [4] and [3] they are as follows. With our selection criteria, the 
typical uncertainty on L(Ho-) is 15% and is dominated by random 
errors. The uncertainty on R* is typically 7%, including a 5% sys- 
tematic uncertainty on the distance modulus. As for the mass M t , 
since it is determined by comparing the location of the star in the 
HR diagram with evolutionary tracks, both systematic and statisti- 
cal uncertainties are important. The uncertainty on the temperature 
is mostly statistical and is about 3%, while that on the luminos- 
ity is 7%, comprising both random errors (1% uncertainty on the 
bolometric correction and 3% on the photometry) and systematic 
effects (5% uncertainty on the distance modulus). When we inter- 
polate through the PMS evolutionary tracks to estimate the mass, 
the uncertainties on T e ff and L* imply an error of 7% on M 4 . In 
summary, the combined statistical uncertainty on M is 17%. 



4 THE NEW SAMPLE OF ACCRETING PMS STAR 
CANDIDATES IN THE LMC 

We have determined the mass accretion rates for 490 PMS star 
candidates in our field #1, 325 in our field #2 and 242 in our 
field #3. Adding the 106 PMS stars in the field around supernova 
SN1987A characterized in Paper I, we end up with a sample of 
about 1000 PMS star candidates in the LMC with mass between 
1 and 5 M Q and age between ~1 and 30 Myr. Our sample rep- 
resents the largest M dataset for low-metallicity stars presented 
so far. The majority of our PMS stars (~90%) have masses be- 
low 2 M , i.e. in the T Tauri star range. The remaining ~10% 
are Herbig Ae/Be stars; at the typical age of our regions (Table[4]l, 
many of them could already approach or be on the main sequence. 
Moreover, they cannot be considered simply as massive counter- 
parts of T Tauri stars because they are characterize d by different 
properties (strong winds, small magnetic fields, etc. : lHerbidll96d ; 
IWaters & Waelkensll 998). Thus, the Ha vs. M calibration used for 
T Tauri stars (Equation[4j might be inappropriate for Herbig Ae/Be 
stars and, hence, the M derived for these objects might be inaccu- 
rate. This does not affect the results presented in this paper because, 
as explained in Sect.[2] we limit our analysis to the range 1-2 M Q . 

HST/WFPC2 photometry, stellar parameters and mass accre- 
tion rates for this sample are available, only in electronic form, in 
Table [2] and [3] respectively. In Table [4] we summarize the number 
of PMS stars selected in each field and their median mass, age and 
mass accretion rate. Using this large sample of PMS stars with M 
measured in a homogeneous way, we now investigate the depen- 



dency of M on stellar mass, age and metallicity, in relation with the 
proposed mechanisms of disk dispersal. 

4.1 On the completeness of our accreting PMS star sample 

In Sect. [2] we estimated that our observations for field #1, #2 
and #3 are photometrically complete in all the filters in the 
mass range 1-2 M Q for 1 to 14 Myr old objects at the dis- 
tance of the LMC and, hence, in this section we limit our anal- 
ysis to this mass range. However, the mass accretion rate of 
PMS stars depends on the stellar ma ss ( Muz erolle et al] 120051; 
iNatta et al]|2006l; ISicilia-Aguilar et al]|2006|) and decreas e s with 
time dMuzerolle et al] l200d ; ISicilia-Aguilar et all I2005L 120061 
l2010l ; lFedele et al]|20ldT Thus, our limit for the measurement of 
M varies within the considered mass and age range. 

Figure [6] shows the M of PMS star candidates in our sample 
as a function of mass (upper panel) and age (lower panel). In each 
panel, the dashed line displays our photometric limit for the detec- 
tion of excess H a emission and, hence, for the measurement of M. 
In order to translate the detection limit in H a luminosity into the de- 
tection limit in M as a function of stellar mass, we considered the 
lowest H a luminosity measured for our PMS stars and used Equ. [3] 
and|4] We assumed the median age of all PMS star candidates in our 
sample (10 Myrs) and calculated the corresponding stellar radii for 
different stellar masses using the 10 Myr isochrone for the metal- 
licity of the LMC. Analogously, we translated the detection limit in 
H tt luminosity into the detection limit in M as a function of stellar 
age by using Equ. [3] and [4] and assuming the lowest H„ luminos- 
ity measured for our PMS star candidates and their median mass 
(1.4 Mo); we then calculated the corresponding stellar radii from 
the 1.4 M Q evolutionary track for the metallicity of the LMC. 

Figure|6]and the limiting magnitudes in TableQ]show that our 
photometric limits allow us to detect only the brightest and, hence, 
youngest accreting late K-type stars at the distance of the LMC, 
while they are sufficient to detect early K or G-type accreting stars 
(Mj > 1.3M Q ) of any age up to the main sequence. In particular, 
the lower panel of Figure [6] shows that our PMS star candidates 
define a bulge of points at ages around 10 Myrs and this bulge is 
very close to the photometric limit; this implies that we might have 
missed part of the PMS objects with very low accretion rates at 
ages around 10 Myr. At younger ages, the mass accretion process 
of PMS stars is enhanced and the H„ luminosity increases; if more 
accreting PMS stars younger than 10 Myr were present, we would 
have them. In summary, our accreting PMS star candidate sample 
in the LMC is complete in the spanned age range for M* > 1.3M Q , 
while we might have missed part of the oldest lower-mass accreting 
stars. 

We now quantify this incompleteness effect by using the 
dataset presented in Paper I as a reference. Indeed, one of the main 
objectives of this work is to compare the accreting properties of 
PMS stars in our LMC fields with the pioneer study presented in 
Paper I, where we found that the M of PMS stars in the region 
around supernova SN1987A are systematically higher than their 
Galactic counterparts (Sect.[6]l. To make sure that this comparison 
is performed in a homogeneous way, we derived a correction to be 
applied to the typical M measured in our LMC fields to match the 
completeness of the Paper I dataset. We proceeded as follows. 

Since the photometric catalog for the SN1987A region (Pa- 
per I) is much deeper than our catalogs, it can be considered as a 
reference catalog with 100% completeness for the measurement of 
M in the considered mass range (1-2 M Q ). We considered the Ha 
luminosity (L(Ha)) distribution of PMS stars in the field around 
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Table 2. Astrometry and observed photometry of the PMS star candidates in the LMC. 



Field 


ID 


RAJ2000 


DECJ2000 


F606W 


F656N 


F814W 






(hh:mm:ss) 


(dd:mm:ss) 








#1 


1 


5:34:47.62 


-69:14:37.89 


22.435±0.028 


20.463±0.109 


21.957±0.028 


#1 


2 


5:34:44.36 


-69:15:13.14 


21.028±0.()12 


19.631±0.051 


20.220±0.()10 


#1 


3 


5:34:43.90 


-69:15:17.41 


21.549±0.015 


20.408±0.092 


21.061±0.015 


#1 


4 


5:34:45.31 


-69:15: 1.14 


19.870±0.008 


18.990±0.027 


19.142±0.007 


#1 


5 


5:34:47.21 


-69:14:41.45 


19.789±0.008 


18.839±0.024 


18.984±0.006 


#2 


1 


5:35:41.65 


-69: 9:10.82 


20.484±0.()02 


18.439±0.017 


19.527±0.002 


#2 


2 


5:35:41.40 


-69: 9: 8.94 


23.139±0.035 


20.382±0.092 


22.338±0.038 


#2 


3 


5:35:41.55 


-69: 9: 8.70 


22.368±0.016 


20.311±0.103 


21.878±0.018 


#2 


4 


5:35:41.69 


-69: 9: 8.25 


22.095±0.009 


20.610+0.111 


21.569±0.012 


#2 


5 


5:35:41.93 


-69: 9: 8.08 


21.438±0.005 


19.892±0.060 


20.738±0.006 


#3-Pop.l 


1 


5:39:32.88 


-69: 9:34.32 


20.324±0.011 


16.878±0.013 


19.452±0.009 


#3-Pop.l 


2 


5:39:36.03 


-69: 9: 8.81 


20.434±0.011 


19.061±0.060 


19.432±0.010 


#3-Pop.l 


3 


5:39:32.47 


-69: 9:40.29 


21.993±0.036 


18.877±0.067 


21.842±0.034 


#3-Pop.l 


4 


5:39:32.94 


-69: 9:36.18 


20.604±0.012 


19.31 1±0.082 


20.331±0.014 


#3-Pop.l 


5 


5:39:35.28 


-69: 9:16.37 


21.993±0.024 


20.400±0.134 


21.354±0.021 



Table 3. Stellar parameters, Ha equivalent width and accretion rate of the PMS star candidates in the LMC. 



Field 


ID 


logT e ff 




Mass 


Age 


EW fl „ 


log M acc 






(K) 




(M Q ) 


(Myr) 


(A) 


(M /yr) 


#1 


1 


3.86±0.02 


0.55±0.09 


1.3 


17.7 


67± 9 


-7.4±0.8 


#1 


2 


3.72±0.01 


1.10±0.04 


2.3 


0.6 


18± 2 


-7.1±0.8 


#1 


3 


3.86±0.01 


0.90±0.05 


1.5 


9.2 


16± 3 


-7.6±0.8 


#1 


4 


3.75±0.00 


1.54±0.03 


3.2 


0.7 


4±0 


-7.2±0.8 


#1 


5 


3.73±0.00 


1.61±0.03 


3.5 


0.3 


3± 


-7.2±0.8 


#2 


I 


3.66±0.00 


1.37±0.02 


1.2 


0.0 


58± 1 


-5.8±0.7 


#2 


2 


3.71±0.02 


0.24±0.12 


1.3 


5.9 


158+16 


-7.2±0.8 


#2 


3 


3.84±0.01 


0.53±0.06 


1.3 


16.4 


73± 9 


-7.4±0.8 


#2 


4 


3.83±0.01 


0.63±0.04 


1.3 


12.0 


31± 6 


-7.6±0.8 


#2 


5 


3.75±0.00 


0.89±0.03 


2.0 


2.7 


32± 3 


-7.2±0.7 


#3-Pop.l 


1 


3.68±0.00 


1.36±0.04 


1.4 


0.1 


292± 5 


-5.3±0.7 


#3-Pop.l 


2 


3.67±0.00 


1.44±0.04 


1.4 


0.0 


17± 2 


-6.4±0.8 


#3-Pop.l 


3 


3.97±0.04 


0.73±0.19 


1.6 


8.5 


245±19 


-7.0±0.8 


#3-Pop.l 


4 


3.91±0.01 


1.21±0.04 


1.8 


6.2 


21± 3 


-7.2±0.8 


#3-Pop.l 


5 


3.80±O.Ol 


0.73±0.08 


1.4 


8.5 


29± 7 


-7.4±0.8 



SN1987A and derived its peak by performing a Gaussian fit (Fig- 
ure[7}. This peak (L(Ha) Re f) represents the limit of this dataset for 
the detection of Ho- excess emission and, hence, for the measure of 
M. 

In Figure [7] we compare the L(Ha) distribution of the 
SN1987A PMS population with the analogous distributions for our 
region #1. This region is very close to SN1987A, has the same age 
(Table [4]l and is part of the same star-forming complex. Thus, he 
difference between the peak of this distribution and L{Ha) Re f rep- 
resent the L(Ha) correction to be applied to field #1 to match the 
completeness of the SN1987A dataset. The L(Ha) correction has 
been translated into M correction by using Equ. [3] and [4] assum- 
ing the median mass and age of the PMS population in region # as 
reported in Tableland using the relative evolutionary model for 
the metallicity of the LMC to calculate the corresponding stellar 
radius. 

The Ho- observations of fields #2 and #3 have roughly the 



same photometric depth as field #1 (Table[T|l and, hence, we assume 
the same M correction. The median M of the PMS star candidates 
selected in each field corrected for completeness are reported in 
Table |4]. 



5 ON THE MASS DEPENDENCY OF THE MASS 
ACCRETION RATE IN THE LMC 

The mass accretion rate of galactic PMS stars appears to be 
roughly proportional to the se cond power of t h e stellar mass 
over a wide range (0.3-3 M ; Nattaetal. 2004; Mohanty et al. 
2005| ; lNatta et al.l2006l : lMuzerolle et alj2005l : rSicilia-Aguilar et al. 
2006), with a large spread of two orders of magnitude at least 
for any interval of M*. Recent observations also show that this 
power law might even be steeper (M oc Af\ A ) in the very low- 
mass and sub-solar mass domain (M t < 0.4M o ), while flatten- 
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Table 4. Number of PMS stars selected in each field and their median mass, age and mass accretion rate. For each parameter the range spanned by each PMS 
population and the standard deviation of the distributions (cr) are also given. 



Field 


N. of PMS stars 


Median Mass 
(Mo) 


Median Age 
(Myr) 


Median log M f 
(M /yr) 


Corrected log M* 
(M Q /yr) 


#1 


490 


1.3 

range [1:5] 
cr=0.5 


13 

range [1: 30] 
tr=6 


-7.5 

range [-8.1: -6.3] 
cr=0.3 


-7.6 


#2 


325 


1.4 
range [1:5] 
<r=0.5 


9 

range [ 1 : 30] 
cr=6 


-7.3 

range [-7.7: -5.3] 
cr=0.5 


-7.4 


#3 - Pop. 1 


96 


1.4 
range [1:4] 
cr=0.4 


8 

range [1: 25] 
cr=6 


-7.1 

range [-7.5: -5.1] 
cr=0.5 


-7.2 


#3 - Pop. 2 


146 


1.5 

range [1:5] 

<T=\ 


6 

range [1:30] 

<T=6 


-7.1 

range [-7.9: -5.3] 
cr=0.5 


-7.2 


SN1987A" 


104 


1.3 

range [1: 2.2] 
cr=0.2 


14 

range [1: 30] 
tr=6 


-7.59 
range [-8.0: -6.4] 
cr=0.3 





"From Paper I. 
^Median mass accretion rate in the mass range 1-2 M Q . 
'"Median mass accretion rate in the mass range 1-2 M Q corrected for completeness effects (Sect. l4~TV 
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Figure 7. Comparison between the Ha luminosity distribution of the PMS 
population in the region around SN1987A (gray histogram, from Paper I) 
and the distribution for our fields #1 (black histogram). The curves show 
the Gaussian fit to the data. 



ing (M oc M\ or e ven flatter) when only the > O.8M do - 
main is considered dFang et alj|2009l ; Isicilia-Aguilar et alj|2010h . 
This seems to support the idea that the accretion and disk evolution 
proce ss might present impo rtant differences for low and high-mass 
stars dHartmann et alj2 006). 

The large spread observed in the M vs. M* dependency partly 
arises from the uncertainty in the determination of M, which can 
be as high as one order of magnitude. However, this relation still 
contains an intrinsic scatter due to the age dependency of M, in- 
dividual disk properties and the variability of the mass accretion 
process, as discussed in Sect. [T] Indeed, mass accretion rates for 
younger stars tend to be higher, as expected from the viscous disk 



evolutionary models {Hartmann et al. 1998). The M versus age de- 
pendency appears to be very steep, with M decreasing by more 
than three order of mag nitude within the first 10 M yr of the stel- 
lar life (see Figure 2 bv lSicilia-Aguilar et al.ll201oT) . This effect is 
likely to increase the spread in M for any interval of , because 
the PMS populations in young clusters and star forming r egions 
in the MW generally present an age spread of a few Myr jjeffriesl 
|20T ]]), and objects progressively younger are expected to form par- 
allel tren ds in the M vs. M» relation with increasing offset in M 
(see, e.g.. lSicilia-Aguilar et al.ll2006l . l201fj|) . 

We now use our sample of PMS stars to investigate the M 
versus M* correlation in the LMC. Although our sample spans a 
wide range of stellar masses (see Table |4j, 90% of the PMS stars 
have masses between 1 and 2 M Q and, thus, we limit our analysis 
to this range, in order to obtain statistically robust results. We find 
that best fit to the M versus M t dependency in the LMC is a linear 
correlation (Figure [6] upper panel). In order to investigate quanti- 
tatively this dependency, disentangling the dependency on age, we 
performed a double linear regression fit, setting M* and age as in- 
dependent variables and M as dependent variable. In a logarithmic 
scale, the functional form of this regression is as follows: 



I02 



M 



a ■ log Age(Myr) + b ■ log 



(5) 



M e yr- [ 

where c is the offset and b and a the coefficients enclosing the 
mass and age dependency, respectively. In this section, we focus 
our attention of the mass dependency of M enclosed in the b pa- 
rameter. The age dependency and the relative a parameter will be 
discussed in the Sect. [6] 

Considering all the PMS stars in the range 1-2 M Q , we find that 
the M of PMS stars in the LMC increases almost linearly with stel- 
lar mass (b = 0.82 ± 0.16, i.e. M ocM+). The solid line in Figure|6] 
(upper panel) displays the results of our linear regression fit; the 
dashed line displays our photometric limit for the detection of ex- 
cess H ff emission and, hence, for the measurement of M (Sect. 14. il l. 
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Figure 6. Mass accretion rate as a function of the stellar mass (upper panel) and age (lower panel) for PMS star candidates in the LMC with masses in the 
range 1-2 M Q . In both panels, the solid line displays the result of the double linear regression fit to the points obtained by setting M* and age as independent 
variables and M as dependent variable, while the dashed line represents our photometric limit for the measurement of M (Sect.[3]and[6). The gray dots are the 
data points in the M* > 1.3M Q mass range, where our PMS star candidate sample is photometrically complete. The dot-dashed line displays the result of the 
double linear regression fit of this sub-set of data. 



In order to avoid the mass regime where our survey is not complete, 
we recomputed the coefficients in Equ.|5]limiting the regression fit 
to the M* > 1.3M range (grey points and dot-dashed line in Fig- 
ure [3J. In this case, we obtain b = 1.02 ± 0.31, confirming that the 
M vs. My, dependency is roughly linear. 



We conclude that in the approximate range 1-2 M Q , the M of 
PMS stars in the LMC increases with the stellar mass as M ocM* 
where b«l, in agreement with the value (b = 0.9 ±0.1) found for 
the PMS population of NGC 346 in the SMC (Paper II). Taking 
into account the uncertainty on the ^-parameter, this dependency is 
slightly flatter than the second power law reported in the same mass 
range for samples of PMS stars in the Galaxy. Using data from the 
IPHAS survey of PM S stars in the galactic HII region IC 1396, 
iBarentsen et al.l d201 lb also found a slope between 1.2 and 1.3 for 



the M vs. M* dependency, i.e. less steep than claimed in previous 
literature. 



6 M EVOLUTION AND METALLICITY DEPENDENCY 

The data collected so far for galactic P MS stars indicate that 
M is a decreasing function of stellar age ( Muz erolle et aI] |200rA 
ISicilia-Aguilar et alj2005Ll2006ll2010l : lFedele et al.l2010l), roughly 



in lin e with the expected evolution of viscous discs i Hartmann et al.l 
Il998l) . alt hough the sp r ead of the data is very large (Sect.[T). More 
precisely. iFedele et alJd2010l) estimated a mass accretion timescale 
(Tacc) of 2.3 Myr for galactic PMS stars in the spectral type range 
K0-M5. Moreover, in the pioneering study of the mass accretion 
process in the MCc presented in Paper I and Paper II, we found 
that, in the field around supernova SN1987A in the LMC and in 
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NGC 346 in the Small Magellanic Cloud (SMC), the M of PMS 
stars are systematically higher than their galactic counterparts. 

Because of the very large uncertainties affecting age and M 
determinations, one must be very careful when doing such a com- 
parison for individual young stars. In particular, isochronal age es- 
timates for PMS stars in the MCs are more uncertain than ages of 
galactic PMS stars, that can be confirmed by other methods (turnoff 
of the main sequence, dynamical evolution of the gas, etc.), and star 
forming regions in the MW are usually s maller and, hence, present 
smaller age dispersion felmegreenl 2000). However, age differences 
based on statistical samples are more reliable. 

In this section, we investigate the M evolu tion in the MCs 
compared to the MW and verify the result by jDe Marchi et al.l 
l20ldl2011ah . We dispose of an enlarged sample of PMS stars in 
the LMC with M measured in a homogeneous way and, hence, we 
adopt a robust statistical approach which is expected to reduce as 
much as possible the effect of uncertainties on age/M. However, 
there are unavoidable caveats that might affect our results and we 
will discuss them in detail in Sect. [7] 

Figure [6] (lower panel) gives a first hint on the time evolution 
of M in the LMC, which, as expected, decreases with age. In or- 
der to investigate quantitatively the M in the LMC we used at first 
the linear regression fit already describe in Sect.[5](Equation[5} and 
shown by the solid line in Figure[6](lower panel). We obtain a ~ 1/3, 
both when considering all the mass range 1-2 M Q and when lim- 
iting the fit to the M t > 1.3M Q regime to avoid incompleteness 
effects. This would imply that the M of PMS stars in the LMC de- 
creases with age significantly more slo wly than observed for Galac - 
tic PMS stars of the same mass (asl ,2; ISicilia-Aguilar et al.l20icil) . 
where the M evolution is marginally consistent, though still slower, 
than the prediction of the viscous disk evolution model (a= 1.4-2. 8; 
lHartmannet"aT]|l998h . This comparison is based on age/M esti- 
mates for single stars and there are at least two main problems af- 
fecting its outcome: 

(i) In Sect. l3.3l we have already pointed out that isochronal ages 
of individual objects are uncertain, especially below a few Myr, 
and only ages and age differences of statistical samples should be 
trusted^ 

(ii) [sicilia-Ag uilar et ail feOld) already pointed out that the dif- 
ference between the observed M evolution and prediction of the 
viscous disk model might be attributable to the fact that solar-type 
objects present a v ariety in the viscosity parameter radial exponent 
(Isella et al. 2009) and long- surviving disks may be biased toward 
certain radial viscosity laws. 

A more robust approach to study the evolution of M in the 
MCs in comparison to our MW consists of deriving the median 
age and M of a sample of star forming regions in both galaxies 
spanning a suitable age range. We have collected from the literature 
17 Galactic star-forming regions and clusters with ages between 1 
and 30 Myr whose stellar population is well characterized in terms 
of M. We considered only Galactic PMS stars with mass in the 
range 1-2 M Q , i.e. the range we are studying in the LMC. This 
sample consists of ~170 objects; their typical age and median M, 
together with the number statistics per age bin and the reference 
papers, are reported in Table[5] 

In Figure [8] we plot the median M as a function of the me- 
dian age for each one of these 18 Galactic PMS populations, for 
the LMC p opulations studied in t his paper and for NGC 346 in 
the SMC foe Marchi et al.l 1201 fal . hereafter Paper II). Note that 
NGC 346 displays two different stellar populations with median 
ages of ~1 (NGC 346-young) and -20 Myr (NGC 346-old), re- 



spectively, and our field #3 in the LMC also presents two spatially 
separated PMS populations (#3-Pop. 1 and #3-Pop. 2) with slightly 
different median ages (6 and 7 Myr, see Sect. 18.4b . Thus, we dis- 
pose of 7 regions in the MCs with ages between 1 and 20 Myr (Ta- 
ble[4j. The errors bars in Figure[8]for the MCs points are computed 
as the typical standard deviation of the M and age distribution of 
each population plus the intrinsic uncertainty on these parameters. 
For the Gala ctic regions we adop t the t ypical error on M and age 
estimated by Sicilia-Aguila r et all d2006h - 

The analysis of Figure[8]reveals that below 3 Myr the only in- 
formation we have at the moment in the MCs comes from NGC 346 
in the SMC. Moreover, in this age range the number statistics of 
Galactic PMS stars with mass between 1 and 2 M Q and measured 
M is very poor (i.e. only ~10% of the Galactic sample). As dis- 
cussed in Sect.[7] below 3 Myr age determination becomes critical. 
Thus, we focus on the age range >3 Myrs and defer the analysis 
of the very young regime (1-2 Myr) to a future paper based on im- 
proved number statistic of 1- 3 Myr old PMS stars in t he core of the 
30 Doradus starburst cluster ( De Mar chi et al J201 lbh . 

After the first 3 Myr, despite the large uncertainties, the av- 
erage M of PMS stars in the MCs is systematically higher than 
that of Galactic PMS stars of the the same age, confirming the 
results by Paper I. Note that the pentagrams in Figure [8] display 
the median M for the two PMS populations in the Orion Nebu- 
lar Cluster (ONC-young and ONC-old, see Table [5]l and the di- 
amond corresponds to Trumpler 37 (Tr 37), for which mass ac- 
cretion rates are phot ometrically derived by Panagia et al. (2011, 
in preparation) and Barents en et al.l feOHl) , respectively, using the 
same method adopted in this paper. These three points follow very 
well the Galactic trend . In particular, the median M measured by 
iBarentsen et al.l J20 1 lh for the PMS p opulation in Tr 37 perfectly 
agree s with the median M measured bv lSicilia-Aguilar et alj J2006I . 
|2010|) for the same population on the basis of U-band excess emis- 
sion and high resolution spectroscopy, which provides accurate 
spectral type and reddening estimate. This demonstrates that the 
difference in the typical M between the MCs and the MW is not 
due to systematic uncertainties of the photometric method we use 
to measure M. 

If we perform a linear fit to these averaged data for ages above 
3 Myrs, we find that in the MCs M decreases with time as M oc 
t~" with a = 1.5 ± 0.2, confirming the average trend in the MW 
(a = 2.6 ± 1.0) and in agreem ent with the prediction of viscous disk 
evolution models (a=1.4-2.8; lHartmann et aljl998l) . 

Being the metallicity of the MCs considerably lower than the 
MW, Figure [8] suggests a tantalizing metallicity dependency of the 
mass accretion process, specifically that M appears to be inversely 
proportional to stellar metallicity. Under this hypothesis, we should 
observe some difference in M even between the two MCs, because 
their metal content is diffe rent (Z=0.007 for LMC and Z=0.002 for 
SMC; lMaederetalJ fl999) and the difference is comparable with 
the difference between the LMC and the MW (Z«0.018). We do 
not observe such a difference when considering the mean M val- 
ues, as we did in our Figure [8] However, Paper II (focused on the 
the SMC cluster NGC 346 and based on individual age and M mea- 
surements) indicates that there is a difference between the two MCs 
and M in the SMC seems to be slightly higher than in the LMC. 
This would fit the idea of M being inversely proportional to the 
metallicity. 

These results suggest that the mass accretion process in the 
MCs might be different from the typical mechanism observed in 
Galactic PMS stars in two main aspects, which are discussed in 
the following sub-sections. In Sect.[7]we will further discuss these 
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Table 5. Median M and number statistics per age bin for galactic PMS populations in the mass range 1-2 M Q . 
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Galactic region 


Typical Age 
(Myr) 


Median log M 

(Ms/yr) 


Number of members with 
1< M/M Q <2 


Ref. 


p-Ophiuchus^ (p-Oph) 


0.5 


-7.95 


4 


Natta et al. (2006) 


Orion Nebula Cluster (ONC-young) 

Orir\n XToKnln finder- (HMP rJHVfc 

w ri u ii iNeouia v^iusLer (win^-okij 


3 

1 z 


-7.60 
-y.uu 


9 

Q 


Panagia et al. 201 1 (in preparation) 


IC348 

T iiniic TTT R nnTTTA 

ijUpus 111 ^i^upiiij 
o"-Orionis (cr-Ori) 
Taurus (Tau) 


2-3 


-7.81 


9 


Dahm (2008) 

ricrLzeji oc niiienurdnu izuuoj 
Rigliaco etal. (2011) 
Calvet et al. (2004) 


Trumpler 37 (Tr37) 
Chamaeleon II (Chall) 


4 


-9.49 


88 


Sicilia-Aauilar et al. (2006. 2010): Barentsen et al. (2011) 
Snezzi et al. (20081 


//-Chamaeleontis (f/-Cha) 

Orion OB lc Association (Ori-OBlc) 

/1-Orionis (/1-Ori) 


5-7 


-8.98 


7 


Javawardhana et al. (2006) 
Calvet et al. (2004) 
Calvet et al. (2004) 


Hydrae (Hya) 


8-10 


-10.52 


4 


Pinzon et al. (2007); Jayawardhana et al. (2006) 


/3-Pictoris (8-Pic) 
NGC 7160 


12 


-10.26 


5 


Pinzon et al. (2007): Sicilia-Aauilar et al. (2010) 
Sicilia-Aguilar et al. (2010) 


Lower Centaurus Crux (LCC) 
Upper Centaurus Lupus (UCL) 


16 


-8.32 


40 


Pinzon et al. (2007) 
Pinzon et al. (2007) 



Tucana-Horologium (Tuc-Hor) 30 -11.03 8 Pinzon et al. (2007): Javawardhana et al. (2006) 



' Note that, while the age of p-Oph is often listed as ~0.5 Myr, Wilking et al. 1 2008) state that the the well characterized, optically visible PMS stars might be 
as old as 2-5 Myr. Thus, in Figure[8]we indicate with an arrow the possible age range spanned by p-Oph members. 
^For the ONC-old population we report the minimum age estimated by Panagia et al. (201 1, in preparation). 



results considering the many caveats of current PMS evolutionary 
models, which might affect age estimates both in the MCs and in 
our MW and hence, our results. We will also review the main mech- 
anisms of disk dispersal (Sect.[8} in order to understand whether the 
lower metallicity of the MCs with respect to the MW plays in favor 
of a typically higher M. 

6.1 Disk lifetime 

For PMS stars in the MW most signs of accretion disappear af- 
ter ~10 Myr. This, together with the observation that the fraction 
of stars with near-IR excess, which originate from the hot dust 
in the inner disk, decrease s at a constant rate jHaisch et al]|200ll ; 
ISicilia-Aguilar et alf2 006). yields an overall disk lifetime of about 
5-6 Myrs. However, a few cases of strong accretors at older ages 
have been reported in the literature. Sicilia-Aguila r et al]j2005h no- 
tice the presence of several Galactic G-type stars, with an apparent 
age older than ~10 Myr, showing higher accretion rates outside the 
typical trend of Galactic PMS stars (see their Figure 12). However, 
they warn that their age estimates might be significantly affected 
by the uncertai nties in the birth lin e for G-type stars (see Sect.[7](. 
More recently, iBaume et al. I J201 lj) also found a handful of stars 
with Hff emission typical of accretors in the 20-30 Myr old cluster 
NGC 6167. 

Moreover, among the 17 Galactic regions considered in Fig- 
ure [8] there is a clear outlyer, represented as an upside-down trian- 
gle: the 16 Myr old Lower Centaurus Crux (LCC)/ Upper Centau- 
rus Lupus (UCL) region, where the typical M is higher than ex- 
pected for its age and is fairly consistent with (though still lower 



than) that measured for the 20 Myr old population of NGC 346 
in the SMC. IPinzon et all j2007l) presented a comparative study of 
the average mass accretion rates in Galactic young clusters and as- 
sociations with ages between 10 and 30 Myr based on photomet- 
ric excesses in the U band; they already noticed that an important 
fraction (~25%) of stars in LCC/UCL presents a relative strong ac- 
cretion. According to these authors, the peculiar case of LCC/UCL 
might indicate that gas reservoirs in disks can exceptionally exist 
much longer than 10 Myrs, allowi ng the formation o f giant planets 
up to the age of ~16 Myr (see also Moor et al. 201 1). However, no 
strong correlation was found between U-band excess emission and 
Ho--excess emission for the stars in the LCC/UCL and some doubts 
arouse about the nature of the U-band excess (R. de la Reza, private 
communication). 

For the sake o f clarity, we also point out that the study by 
IPinzon et all J2007t) finds 8 accretors in Tuc-Hor and 5 in /3-Pic in 
the mass range 1-2 M Q , wh ile previous studies of the PMS popu- 
lation of these associations d Javawardhana et aU l2006) find only 2 
accretors in /?-Pic and none in Tuc-Hor. The authors do not discuss 
this discrepancy. 

In the MCs, we measure an average M of 5-10~ 8 M o /yr for 
10-12 Myr old stars, indicating that disks in the MCs might be 
long-lived with respect to the MW. Observational measurements 
of the disk lifetime in the LMC and, in general, low-metallicity 
environments are at the moment very poor and controversial. Us- 
ing data from the Spitzer SAGE Survey in the MCs, it was found 
that the fraction of stars with disks in the SMC is higher than in 
the LMC (D. Lennon, private communication). Being the SMC 
more metal poor than the LMC, this would be consistent with a 
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Figure 8. Median mass accretion rate as a function of the median age for PMS populations in the LMC (filled circles), NGC 346 in the SMC (filled triangles) 
and 17 young clusters and star forming regions in the MW (open squares), as indicated in the labels. T he upside-down triangle represents the measurement 
for the LCC/ UCL region, while the diamond and the pentagrams are the photometric measurements by Barentse n et alj 4201 ll) and Panagia et al. (2011, in 
preparation), respectively. The arrows indicate the age uncertainty for the p-Oph members and the ONC-old population (see Table|5J. The typical errors on M 
and age are displayed in the top-left and bottom-left for the MC and MW regions, respectively. 



scenario where the disk lasts longer in a metal poorer environ- 
ment. On the other hand, in the extreme outer Galaxy (E OG), where 
the me tallicity is one-tenth of the solar neighborhood, lYasui et ail 
(2009) found a fraction of stars with circumstellar disks signifi- 
cantly lower than in the solar neighborhood, suggesting that most 
stars in low-metallicity environments experience disk dissipation at 
earlier stages (< 1 Myr). The mism atch between these two results is 
rather awkward. Yas ui et alj d2009h might have underestimate the 
disk fraction in the EOG because their measurements are based 
only on near-IR excess, w hich can be very s mall or inexistent in 
flattened disks. Moreover, lYasui et alj (2009) warn that a mecha- 
nism specific to the EOG environment, i.e. effective far-ultraviolet 
phot oevaporation, m ight contribute to the rapid disk dispersal (see, 
e.g.. lHavwoodll20?J8l) . The mechanisms of disk dispersal and their 
dependency on metallicity are further discussed in Sect. [8lto |8.4.2l 



6.2 Disk mass 

In order to keep an average M of 5-10~ 8 M Q /yr for ~10 Myr, the 
disks we are observing in the MCs should be very massive, i.e. 
about 0.5 M G . Moreover, it is normally assumed that accretion 
rates are much higher in the protostellar phase (<1 Myr) than at 
5-10 Myr. Thus, if it is really true that PMS stars in the LMC have 



average M of 5-10~ 8 M Q /yr at 10 Myr, one would expect them to 
have even higher M at earlier ages and, hence, even higher disk 
mass. 

Disk properties and initial conditions may be different in low- 
metallicity environments and disk dispersal mechanisms other than 
accretion might c ontribute to diss ipating the disk mass. Moeover, a 
recent study bv lZhu et alj (I2OIO) suggests that, if the initial cloud 
core is moderately rotating, the protostar ends up with a more mas- 
sive disk and most of the central star's mass is built up during FU 
Ori-like outburst events of accretion. However, it is still hard to ex- 
plain how such massive disks are produced, because oservations in 
the MW suggest that stars accrete most of their mass in the embed- 
ded phase, and only a few percentage once they reach the T Tauri 
phase. 

The high M measured in the MCs would imply that the overall 
accretion process during the PMS phase will add to the central ob- 
ject a mass comparable to that of the central object itself and, hence, 
the protostars itself should have as much mass in the central object 
as in the disk/envelope. This appears unlikely, because the Initial 
Mass Function (IMF) in the MCs appears remarkably sim ilar to the 
that of the MW (e.g. jGarmanvlll998l : lBastian et alj20l"ch . 

However, there is an important point to be considered when 
doing such speculations. The mass accretion process of stars in 
the PMS phase is a highly variable phenomenon and this vari- 
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ability is reflected in the diagnostics used to measure M ( Ha ex 



cess emission, veilin g, mid-IR Bry and Pa/3 em ission, etc.; Herbst 
19861; iHartigan et jdll99ltlNpven et al]l2009l ; lFedele et alj|20ld ; 
Petrov et alj201 luFaesi et alj201 11) . At any given time, only a cer- 
tain fraction of PMS stars of the given PMS population shows sig- 
natures of mass ac cretion above the de tection limit of the given sur- 
vey. The study by iFedele et ail d201Ci) , based on Ho- excess emis- 
sion measurements, shows that the fraction of PMS stars with M 
above their detection limit in Galactic star forming regions de- 
creases with age. At the typical age of our PMS star candidates 
(~10 Myr), the fraction of PMS stars with M above their detection 
limit is expected to be of the order of 5% (see Figure 3 and 4 by 
Fedeleetal j201Ch or, in other words, their M is above the detection 
limit only for 1/20 of the time. 

Since the mass accretion process depends on metallicity (see 
Sect. [HJ, this estimat e might not be accura te for low-metallicity 
PMS stars. Indeed, |Pe Marchi et al.l j201ld) found that the frac- 
tion of PMS stars above their M detection limit at ages younger 
than ~ 10 Myr in the SMC cluster NGC 346 is of the order of 30%. 
In other words, these PMS stars are above the M detection limit 
only for 1/3 of the time. The metallicity of the LMC is in between 
the SMC and the MW; we estimate that our PMS star candidates 
in the LMC, with an average M of 5-10~ 8 M G /yr, a typical age of 
~10 Myr, accrete a total mass of ~0.02 up to ~0.2 M Q depend- 
ing on whether we assume that their are above our M detection 
limit (dashed line in Figure |6]( 1/20 or 1/3 of the time. This im- 
plies that the disk mass is between 10 and 20% of the stellar mass, 
because our PMS star candidates have mass in the range 1-2 M B , 
which is plausib le dWhite & Hillenbrandl2004l ; lAlcala et al.ll2008l ; 
iKim et al.ll2009h . 



7 CAVEATS ON OUR AGE AND M ESTIMATES 



As mentioned in Section [331 the uncertainty on stellar masses and 
ages derived from PMS evolutionary models could have a signifi- 
cant impact on our results. In particular, a few caveats on isochronal 
ages should be taken into ac c ount. 

Recently iBaraffe et al] |2009) showed that, when vigorous 
episodes of mass accretion at early stages of the stellar life 
(< 1 Myr) are taken into account in the calculation of PMS evolu- 
tionary tracks, protostars may have much smaller radii than found 
in previous treatments. Such small radii would have the effect of 
making some young stars appear fainter and thus much older in the 
HR diagram, perhaps as much as 10 Myr, than they really are. This 
means that the concept of stellar birth-line has no valid support for 
very low-mass stars because it depends on the accretion history. 
As a consequence, ages below a few Myr are highly uncertain and 
absolute ages are hard to determine, though relative ages are still 
reliable. This problem has a strong impact on M evolution stud- 
ies, because it implies that a strongly accreting object might show 
a position on the HR diagram consistent with a 10 Myr age, even 
though its true age is much you nger. 

However, the prediction by Bar affe et al.l {2009) appli es to the 
limit o f low-temperature ("cold") accretion, and Hartmann et al. 
J20T 1) argue that very rapid disk accretion is unlikely to be cold, 
for several reasons: i) the FU Ori objects, e.g. the best-studied PMS 
disks with rapid accretion outbursts, have spectral energy distribu- 
tions consistent with large, not small, radii; ii) theoretical models 
indicate that at high M, protostellar disks become internally hot 
and geometrically thick, making it much more likely that hot ma- 
terial is added to the star; iii) the luminosity of the accretion disk 



is likely to irradiate the central star strongly, heating up the outer 
layers and potentially expanding them, iv) observed HR diagram 
positions of most young stars are inconsistent with the rapid cold 
accretion models. 

The IBaraffe et alj (2009) prediction implies that we might 
have overestimated the age of our PMS star candidates. This prob- 
lem affects age estimates both in the MCs and in our MW; however, 
under the hypothesis that the typical M are higher in the MCs, it is 
expected to be more significant for PMS stars in the MCs. In other 
words, if we accept a typically higher M in the MCs, then we have 
to accept that our age estimates might be inconsistently derived, 
because isochronal ages are not robust for strongly accreting ob- 
jects and we might have overestimated the age of PMS stars in the 
MCs. This problem can be properly addressed only when improved 
isochrones for low-metallicity PMS stars will be available. How- 
ever, this issue does not affect considerably our conclusion based 
on Figure [8] because the IBaraffe et alj (2009) effect is expected to 
be more significant for very young objects (<1 My) of spectral type 
M or later, ranging from a few Jupiter masses to a few tenths of a 
solar mass, while our sample consists essentially of G and K type 
stars with ages older than ~5 My (see Table[4]l. 

Hartmann (2003) discussed the uncertainties in the birth line 
of G-type stars (T f // >5400 K), which results in a overestimate of 
their ages and dominate the apparent age spread. We do observe 
that G-type star in our sample are systematically older than K-type 
stars by a few Myr. This age difference is partly explained by the 
fact that our photometric limits (TableQJ allow us to detect only the 
brightest and, hence, youngest late K-type stars at the distance of 
the LMC, while they are sufficient to detect G-type stars of any age 
and spectral sub-class up to the main sequence (see also Figure [6] 
upper panel). However, we can not exclude some residual effect due 
to the uncertainty in the birth line of G-type star. This would point 
toward an overestimation of the typical age of PMS populations in 
the MCs. 

Finally, in Sect. 13.31 we pointed out that the regions we are 
studying are very large, about 100 pc in diameter, and we expect 
an age dispersions of a few to 10 Myr within each region (Tabled. 
This age dispersion results from the presence of multiple stellar 
populations in the same region, which mi ght be the resu lt of se- 
quential and/or triggered star formation JElmegreenll 20001) . A typi- 
cal example in our MW is the Cep OB2 bubble, which is ~50 pc in 
diameter and includes the 12 Myr old cluster NGC 7160, the 4 Myr 
old Tr 37, the 1 Myr old Tr 37 globule, the 2-3 Myr old Cep B, and 
several other star formation episodes with ages ranging 1-3 Myr 
dPatel et al.lll998l ; ISicilia-AguiIar et ai1l2005l) . 

Similarly, the three LMC regions analyzed in this work present 
a significant age spread. Figure [5] shows the HR diagram for the 
PMS population of these regions and indicates that they have a typ- 
ical age of ~13 Myr, ~9 Myr and 6-8 Myr for region #1, #2 and 
#3, respectively; however, a small fraction of the PMS population 
in each region is as young as 1 Myr. Since it is normally assumed 
that accretion rates are much high er at ages around 1 Myr than at 
5-10 Myr (see, e.g., Figure 2 by ISicilia-Aguilar et alj|201 0). the 
median M we measure could be driven by the younger populations 
present in each region. This problem affect also the age and M es- 
timates for the MW regions, but it is expected to be less significant 
because these regions are much smaller. 

On the other hand, our conclusion in Sect. [6] are driven by 
Figure [8] using the median age and M of each region in the MCs 
and the MW. This statistical approach is expected to smoothen the 
age uncertainty due to strong accretors and/or G-type stars and, in- 
deed, the median isochronal ages estimated for the LMC regions 
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are supported by several arguments and independent age estimates 
(Sect. 13.31 ). However, considering all the caveats listed above, it 
could be that we are slightly overestimating the age of the MCs re- 
gions. Under the most conservative hypothesis that the MC regions 
are a few Myr younger than our estimates, we still find a systematic 
and interesting differences between the MCs and the MWs, as the 
accretion rates in the MCs are about an order of magnitude higher 
than in Galactic regions of similar age. As we will see in Sect. [8] 
this difference can be justified on the basis of the different metallic- 
ity between the MCs and the MW. 



8 MECHANISMS OF DISK DISPERSAL : WHY 

LOW-METALLICITY DISKS MAY HAVE HIGHER 
ACCRETION RATES? 

Our pioneer study of the mass accretion process in the LMC sug- 
gests that disks around metal-poor stars accrete at very high rates 
with respect to galactic PMS stars with roughly solar metallicity. 
The many caveats of age/M measurements do not allow us to quan- 
tify such a difference with high accuracy. However, all our argu- 
ments support a difference of at least an order of magnitude in M. 
Models of disk evolution provide several explanations both for and 
against this observational evidence. In the next sections we review 
the interplay among the main mechanisms of this dispersal in low- 
metallicity environments and in relation to our finding in the LMC. 



8.1 Viscous accretion 

Viscous accretion is the mechanism by which mass is accreted onto 
the central star because the transport of angular momentum out- 
ward allows disk material to flow radially inward ( Holl enbach et al.l 
200G). In general, a lower metallicity implies a lower disk opacity, 
lower disk temperature, l ower viscosity, and thus longer viscous 
time (see, e.g., review bv lDurisen et al.l 2007). Therefore, it is ex- 
pected that low-metallicity stars will undergo significant accretion 
up to older ages with respect to higher metallicity stars, in general 
agreement with our measurements in the LMC. However, a few 
warnings must be taken into account: 



(i) According to dHartmann etai] l2006h , in disk around low- 
mass stars the accretion/viscosity is mainly driven by magneto- 
rotational instability (MRI), which may be more effi cient in low 
metall icity stars. Indeed, this is the argument used by I Yasui et al.l 
(2009) to explain the fact that in the EOG, where the metallic- 
ity is one-tenth of the solar neighborhood, the fraction of stars 
with circumstellar disk is significantly lower, implying a faster disk 
dispersal. This consideration would go against he long-lived low- 
metallicity disks hypothesis. Nevertheless, there is no observational 
evidence that accretion/viscosity driven by MRI is so efficient. For 
example, evolved disks with large grains, which are expected to be- 
have like low-metallicity disks, accrete less as they age, while, un- 
der the assumption that the MRI is more effic ient at low-metallicity , 
we would expect the opposite (see Sect. 5 bv lHartmann et alj2006l 
and references therein). 

(ii) In order to keep an average accretion rate of 10~ 7 M Q /yr 
going on for 5-10 Myr, the disks we are observing in the LMC 
should be extremely massive, about 0.5-1 M Q . An additional ef- 
fect to take into account in massive disks is gravitational insta- 
bility (GI). It has been suggested that lower opacity favor GI 
dCossins et a l. 2009) and that low m etallicity disks may be more 
unstable against GI dCai et alj r2006). There is an ongoing debate 



about whether GI makes ve ry massive disks fragment dBossll 19971 
120031 : 1 Cai ct al ]|2006l 2009), or simply favor accretion and mass 
transport dHartmann et alj2006h . In both cases, the GI global effect 
would point against the long-lived low-metallicity disk hypothesis, 
even though it could explain the large M we measure. 



8.2 Grain growth 

According to the analysis by iDullemond & Dominikl d2005l) . the 
growth of dust particles in the disk depends on the disk's dust to gas 
ratio. Low-metallicity disks have a lower dust content and, hence, 
the process is slower and the grains remain small. Grain growth is 
an important contributing factor to disk dissipation by planet for- 
mation and photo-evaporation, which we will analyze later on. In 
addition, if the disk is less optically thick, the convective turbu- 
lence is reduced, making the grain coagulation/fragmentation pro- 
cess less efficient. In other words, low-metallicity disks may reach 
a quiet steady state characterized by small grains that cannot co- 
agulate and, hence, they last longer, because dissipation by planet 
formation and/or photo-evaporation is less efficient. This scenario 
is in agreement with our measurements in the LMC, but again there 
are a few warnings to take into account: 

(i) Convection is not the only source of turbulence. As seen in 
Sect. 18. l| MRI driven turbulence will be present in the magneti- 
cally active parts of the disk. However, the high M we measure 
at older ages suggests large disk mass to sustain such long accre- 
tion lifetimes and, hence, a high gas column de nsity. This kind of 
disks would present m agnetically "dead zones" (Tur ner et al J 20071 : 
iTurner & Saiicl [2008) with lower turbulence. These dead zones 
might suffer from other types o f instability (e.g. streaming insta- 
bility, Johansen & Youdin 2007) and they may be good locations 



for planetesimal formation. However, as seen in Sect. |6land [8Tl 
such large disk masses are difficult to explain. 

(ii) The lower dust content of low-metallicity disk implies 
that the gas reservoir of the disk is less shielded against photo- 
evaporation and it might depleted faster. Disk photo-evaporation ef- 
fects and their dependency on metallicity are discussed in Sect. |8.4| 



8.3 Planet formation 



Ijohansen etafl d2009ri presented three dimensional numerical sim- 
ulations of particle clumping and planetesimal formation in proto- 
planetary disks with varying amounts of solid material. These sim- 
ulations showed that the formation of planetesimals is strongly de- 
pendent on the dust to gas ratio. Low-metallicity disks would not 
be able to produce planets at the beginning, but they may do it later 
on once other mechanisms, like internal photo-evaporation, have 
reduced the amount of gas. Thus, the correlation between host star 
metallicity and exoplanets may reflect the early stages of planet 
formation and low-metallicity disks can be particle enriched dur- 
ing the gas dispersal phase, leading to a late burst of planetes imal 
formation. The metallicities studied by Ijohansen et al.l d2009l) are 
solar or super-solar, not as low as in the LMC. However, if low- 
metallicity disks would produce planets only in late phases, they 
might be able to keep an average M higher than higher metallicity 
disk, because early proto-planet growt h might reduce the accretion 
to the central star dMerin et alj|2010l . and references therein). In 
a very low-metallicity environment such as the LMC, this effect 
might be enhanced and would explain the large M we measure. 
For the sake of completeness, we defer the reader to the fol- 
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lowing works pointing out other open issues and possible bias on 
the planet-metallicity correlation: 



• lYasui et al, found that most stars in low-metallicity en- 
vironments experience disk dissipation at earlier stage (<1 Myr) 
than in the solar neighborhood and argue that this could be one of 
the major reasons for the strong planet-metallicity correlation (see 
Sec. [I), because the shorter disk lifetime reduces the time available 
for planet formation. This result is in disagreement with our mea- 
surements for low-metallicity PMS stars in the LMC, whose accre- 
tion process appears to last longer than for galactic PMS s tars. In 
Sect. [6] we have seen that the results by lYasui et al.l |2009) might 
suffer incompleteness effects and the rapid disk dispersal in the 
EOG might be due to effective far-ult raviolet photo-evaporation 
rather than planet formation. Moreover, Ercolan o & Clarke] i201Ch 
explored the metallicity dependence of a dispersal mechanism 
based on planet formation in the core accretion scenario dlda & Linl 
1 20041) and found opposite results, i.e. planet formation in low- 
metallicity environments would extend disk lifetime, because the 
lower metal content implies a smaller amount of dust and, hence, 
a slower planet formation process. These authors explain the ob- 
servations in the EOG in the context of a disc dispersal mecha- 
nism based on X-ray photo-evaporation from the central star, which 
operates faster in low metallicity environments (see discussion in 
Sect.[ia. 

• According to iBouchv et al. (2009), planets formed in lower- 
metallicity disks would accrete less gas and would not suffer so 
much migration; thus, they may end up as low-mass planets very 
hard to detect. 



8.4 On the disk dispersal through photo-evaporation 

Photo-evaporation of the surface layer of the gas in the disk can be 
triggered by strong radiation fields generated by nearby OB stars 
(external photo-evaporation) and/or the radi ation field of the centra l 
PMS star itself (internal photo-evaporation lHollenbach et al. 2000). 
These two mechanisms operate very differently and have very dif- 
ferent responses to changes in metallicity. Thus, we analyze them 
separately. 



8.4.1 Internal photo-evaporation 

Internal photo-evaporation is the mechanism by which radiation 
from the central source heats the gas in the disk, enabling it to 
climb up the gravitational potential well and escape. Theoreti- 
cal studies of the erosion of circumstellar disks around solar- 
type stars via internal phot o-evaporation suggest typical mass 
loss rates of ~ 10~ 10 M o /vr dAlexander et al.ll2006 3 lbh for photo- 
evaporation caused by extreme ultra-violet (EUV) r adiation and 
~10~ 8 M o /yr for photo-evaporation caused by X-rays dOwen et"al] 
l20ld) . Thus, X-rays appear to be much more efficient at drivin. 



internal photo-evaporation (see discussion in Ercolano et alJl2009T ; 
lErcolano & Owenll20ld) . lErcolano & Clarke! J2010h explored the 
metallicity (Z) dependence of X-ray mass loss rates and found that 
X-ray wind rates are higher in lower-Z environments and, con- 
sequently, disk lifetimes are shorter. This occurs because heavy 
elements represent the main source of opacity for X-rays; the 
low metallicity implies a reduction of the X-ray opacity and, 
hence, a larger column c an be heated by X-rays (see Figure 2 of 
lErcolano & Clarkell2010h driving a denser wind from deeper down 



in the dislfl We find that metal-poor stars in the LMC accrete at 
very high rates with respect to galactic PMS stars of the same 
age, whic h can be explained in th e context of the models pre- 
sented by lErcolano & Clarke! d2010h . These authors provides pre- 
liminary evidence for shorter disk lifetimes at lower metallicities 
when X-ray driven photo-evaporation is considered as the main 
disk dispersal mechanism. To a first approximation a disc will be 
destroyed by photo-evaporation when M=M W (Z), where M is the 
accretion rate through the disk and M W (Z) is the X-ray wind rate. 
lErcolano & Clarke! BoTch show that M w oc Z" 77 . This depen- 
dency is steep enough to compensate for the about one order of 
magnitude difference in M between the LMC and the MW, because 
the LMC has about half solar metallicity. Moreover, mass accre- 
tion rates al lower metallicity are expected to be higher because 
the disk viscosity parameter (a acc ) is higher at lower Z. In the MRI 
scenario this occurs because, thanks to the lower X-ray opacity, a 
larger mass fraction of the disk is ionized and therefore MRI active. 

Thus, a disk dispersal scenario based on internal X-ray photo- 
evaporation is consistent with the idea of higher M at lower metal- 
licity. 



8.4.2 External photo-evaporation 

External photo-evaporation, caused by nearby massive stars, works 
primarily via FUV heating of the gas in the outer disk and relies 
on photoelectric heating from dust grains. The most striking exam- 
ple of external photo-ev a poration is seen in the proplyds in Orion 
dBallv et al.ll 1 998l l2000h ■ iLinskv et al] d2007t) also argued that the 
ultra-violet (UV) emission and winds of nearby hot stars play an 
important role in t he development of evaporating gaseous globules 
dHesteretalll 19961) into protostars. External photo-evaporation is 
difficult to quantify due to a number of uncertainties, including pro- 
jection effects, and there still is no clear indication of its global 
effect on PMS populations. 

To study the effect of external photo-evaporation on the PMS 
populations in our fields, we analyze the spatial distribution of the 
PMS stars and their level of Ha luminosity with respect to the dis- 
tance from the hot OB stars in the field, following the same ap- 
proach as in Paper I. To identify OB stars in our fields, we used 
the catalog of UBVI stellar photometry of bright stars across the 
ce ntral 64 deg 2 a r ea o f the LMC and the relative extinction map 
by IZaritsky et al.l d2004l) . We compare the dereddened photome- 
try of these brig ht stars with the tabulation of absolute magnitudes 
of OB stars b y IWegnerl hood) , assuming a dista nce modulus of 
51.4+1.2 kpc dPanagialll99ll ; IPanagia et alj|l99ll) . Only hot stars 
with spectral type earlier than B 1 were considered, because cooler 
stars do not contribute significantly to the ionization (see Table 1 
in Paper I). The center of mass of the hot stars in each field was 
computed using as weight their bolometric luminosity. 

Figure [9] shows the spatial distribution of the PMS stars in 
our LMC fields and the position of the ionizing stars. The center 
of mass of the ionizing flux is represented in each region by the 
large square. Symbols of different colors indicate PMS stars with 
different level of Ha luminosity, as indicated in the caption, and 
the lines display the contours of the Ha nebular emission retrieved 



3 Note that photo-evaporation driven by EUV or far ultra-violet (FUV) ra- 
diation will not give this result. In the case of FUV, the lack of metals and, 
hence, dust would actua lly hinder hea ting although the reduction of opac- 
ity is still present (Gorti & Hol lenbach 2009), while in the case of EUV 
radiation the main source of opacity is Hydrogen rather than metals. 
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Figure 9. Spatial distribution of the PMS stars in our LMC fields: red 
squares correspond to PMS stars with L(Ha-)> 0.03L o , white squares in- 
dicate objects with L(Har)< 0.01L o , and green diamonds are used for inter- 
mediate values. The blue stars indicate the position of the ionizing stars in 
each field and the large yellow square indicates the relative center of mass. 
The contours of the Ha nebular emission are also drawn (black lines). The 
two big circles in the lower panel indicate the two possible PMS popula- 
tions in field #3 and the two yellow squares show the position of the center 
of mass of the ionizing stars for these two populations. 



from our WFPC3/f656n images. Young stars tend to concentrate in 
region of high dust/gas density and indeed, most of the PMS stars 
in our sample follow the lanes of the nebular emission. We note, 
in particular, that our field #3 presents two concentrations of PMS 
stars corresponding to two well defined over-densities in the neb- 
ular emission. Moreover, one of the two groups is also associated 
to a clear cluster of OB stars. We do not observe a significant age 
difference between these two groups (Tabled, also because of the 
limitations in the age determination (Sect. |3.3t . However, the pecu- 
liar spatial distribution suggests that two populations of PMS stars 
separated by about 0.035 deg (i.e. ~30 pc at the distance of the 
LMC) might be present in this region, as indicated by the circles 
in Figure [9] Thus, in our analysis we deal with these two groups 
separately and refer to them as field #3-Pop.l and field #3-Pop.2. 

Figure [9] shows a variety of behaviors in the spatial distribu- 
tion and typical Ha luminosity of the PMS stars with respect to the 
center of mass of the OB stars for the three LMC fields. To show 
this difference quantitatively, we plot in Figure [10] the frequency 
of PMS stars and their typical Ha luminosity as a function of the 
distance to the closest ionizing OB stars for 5 distance bins. The 
dashed lines display the linear fit to the points, which gives and in- 
dication of the global trend of the data. The same plots have been 
presented in Paper I for the field around the supernova SN1987A 
(see their Figure 13 and 14). The result of this analysis for the LMC 
fields can be summarized as follows. 

For the field around supernova SN1987A, we found in Paper I 
that PMS stars in the vicinity of OB stars are both less numerous 
and fainter in Ho- emission than farther away. This suggests that 
their circumstellar disks have been considerably eroded, and made 
less efficient, by enhanced photo-evaporation caused by the ioniz- 
ing radiation of the massive stars integrated over their ~2 Myrs life- 
time. As shown in Figure[l0] we find a similar behavior in our field 
#3-Pop. 1, where the number of PMS star candidates in the vicinity 
of OB stars is lower than farther away. However, the Ha luminos- 
ity of this PMS population does not show any specific trend with 
respect to the position of OB stars. 

The inspection of Figure [10] for the PMS population of our 
field #1 and field #3-Pop.2 shows a different trend. In these two 
fields there are more PMS stars near the ionizing OB stars, while 
their typical Ha luminosity is roughly constant. These correlations 
are consistent with a scenario where photo-evaporation due to ion- 
izing radiation of nearby OB stars does not play a significant role 
for disk dissipation. 

Moreover, the fact that the distribution of low-mass stars fol- 
lows the distribution of massive stars suggests that i) either they 
belong to the same generation or ii) they formed during separate 
events of star formations sharing the same center of mass. Regions 
#1 and #3-Pop.2 are 11 and 6 Myr old, respectively (Tabled; al- 
though we could have slightly overestimated the age of these re- 
gions (Sect. [TJ, both of them are too old to assume that the low- 
mass PMS stars and the OB stars belong to the same generation. 
Thus, our data support hypothesis ii), i.e. two events of star forma- 
tion separated by a few Myr might have occurred in these regions. 
Indication of multiple events of star formation wit hin the same clus 
ter have been already reported in the literatu r e jPatel et al. 



Beccari et alj|2010l ; 
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The inspection of Figure|9]and[T0]for field #2 does not reveal 
any clear trend. 

Thus, the LMC stellar populations show a variety of behav- 
iors with respect to disk photo-evaporation due to external radia- 
tion fields. Although this result might appear puzzling, there is a 
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number of considerations which help illustrating the complexity of 
the problem and understanding the "no trend" detection in field #2. 
Indeed, even in our Galaxy, it is very common to see disk photo- 
evapora tion effects in indiv idual case (see, e.g., the case of Trum- 
pler 37; lMercer et al . 2009), while finding a "global" trend is much 
harder for several reasons. 

The first reason in purely statistical. Even in a large, massive 
cluster there will be few objects close enough to the massive stars 
(~1 pc) to significantly suffer photo-evaporation effects. In addi- 
tion, we can only measure projected distances and this implies that 
lots of contaminants, that are not as close to the OB stars as we 
might think, are included in the analysis. 

Secondly, low-mass PMS stars need to spend a few Myr within 
1-2 pc distance from a OB star to significantly suffer disk photo- 
evaporation effects, though the timescale of this process is still very 
uncertain. With a typical di spersion velocity of 1-2 km/s (see, e.g., 
iKraus & Hil lenbrand 2008, and reference therein), many low-mass 
stars move 1-2 pc per Myr or more if they are ejected and, hence, 
even if they were originally close to a massive star they may end up 
somewhere else. 

Furthermore, the structure/orientation of the given young clus- 
ter can be such that it happens to have another young population 
in the close foreground/background. This second population is not 
physically connected with the cluster under study, but it is indis- 
tinguishable because we can only measure projected distances. In 
these cases, any correlation with the position of OB stars would be 
even harder to detect. 

Finally, triggered star formation may mimic the effect of ex- 
ternal photo-evaporation. If the winds of a massive stars triggers 
more star formation in the outskirts of the cluster , as seen in the 
galactic cluster Tr 37 ( Sicilia- Aguilar et al. 200^), the disk frac- 
tion is higher at larger distances from the massive star. This occur 
just because the stellar population there is younger and is not a disk 
photo-evaporation effect. 

Because of these overall uncertainties, we consider it prema- 
ture to discuss the metallicity dependency of disc dispersal through 
external photo-evaporation. 



9 CONCLUSIONS 

We presented a multi-wavelength study of three star forming re- 
gions in the LMC spanning the age range 1-14 Myr. We aim at 
investigating the mass accretion process in an environment with 
a metallicity significantly lower than our MW. We identify about 
1000 PMS star candidates actively undergoing mass accretion in 
these regions on the basis of their Ha emission and estimate their 
mass, age and mass accretion rate. Our measurements represent 
the largest M dataset for low-metallicity stars presented so far. The 
conclusions of our study are as follows: 

• In the mass range 1-2 M Q , the M of PMS stars in the LMC in- 
creases with stellar mass as M ccM* with foal, i.e. slightly slower 
than the second power low previously reported for galactic PMS 
stars in the same mass regime, and in agreement with the recent re- 
sults by B arentsen et al.l J20 111) for the galactic HII region IC 1396 
(£=1.2-1.3). 

• We find that the typical M of PMS stars in the LMC is higher 
than for galactic PMS stars of the same mass, independently of their 
age. Considering the caveats of isochronal age and M estimates, the 
typical difference in M between the MCs and our MW appear to be 
about an order of magnitude; 
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Figure 10. Fraction of PMS stars with respect to all stars (left panels) and 
relative median Ha luminosity (right panels) as a function of the distance 
from the closest ionizing star for our fields #1, #2 and the two possible 
PMS populations in field #3 in the LMC. The dashed line in each panel is 
the linear fit to the points. 



• Currently available models of disk evolution/dispersal support 
the hypothesis that the higher M measured in the LMC might be a 
consequence of its lower metallicity with respect to our MW; 

• The sample of PMS populations studied in the LMC shows a 
variety of behaviors with respect to disk photo-evaporation due to 
external radiation fields. In the region around supernova SN 1987 A 
(Paper I) we found clear evidence that circumstellar disks have been 
eroded by photo-evaporation caused by nearby massive stars. In the 
three regions analyzed in this work, there is no clear evidence of 
such an effect. However, the analysis of the spatial distribution of 
the PMS populations in region #1 and #3-Pop. 2 with respect to 
the massive OB stars in the field revealed signs of sequential star 
formation. 

In order to clarify these issues, improved evolutionary mod- 
els for low-metallicity PMS stars and modeling/calculations of disk 
dispersal at low-metallicity are needed, as well as observations of 
slightly older (e.g. 30-50 Myr) low-metallicity regions than those 
presented in this paper. 



ACKNOWLEDGMENTS 

We are grateful to an anonymous referee whose suggestions have 
helped us to improve the presentation of this work. We thank A. 
Natta, L. Testi and G. Beccari, for the many useful discussions and 
suggestions, and M. Romaniello for making available his tool for 



20 L. Spezzi et al. 



calculating stellar masses and ages. This publication makes exten- 
sive use of data products from the HST archive at the Canadian As- 
tronomy Data Centre (CADC), operated by the National Research 
Council of Canada with the support of the Canadian Space Agency. 



REFERENCES 

Adams, R, Hollenbach, D., Laughlin, G., & Gorti, U. 2004, ApJ, 
611,360 

Alcala, J.M., Spezzi, L., Chapman, N., et al. 2008, ApJ, 676, 427 
Alexander, R.D., Clarke, C.J., & Pringle, J.E. 2006, MNRAS, 369, 
216 

Alexander, R.D., Clarke, C.J., & Pringle, J.E. 2006, MNRAS, 369, 
229 

Bally, J., Sutherland, R. S., Devine, D., & Johnstone, D. 1998, AJ, 
116, 293 

Bally, J., O'Dell, C. R„ & McCaughrean, M. J. 2000, AJ, 119, 
2919 

Baraffe, I., Chabrier, G., & Gallardo, J. 2009, ApJ, 702, 27 
Baggett, S., et al. 2002, in HST WFPC2 Data Handbook, v. 4.0, 

ed. B. Mobasher, Baltimore, STScI 
Barentsen, G., et al. 2011, MNRAS, 415, 103 
Bastian, N., Covey, K. R., & Meyer, M. R. 2010, Ann. Rev. A&A, 

48, 339 

Baume, G., Carraro, G., Comeron, E, & de Ella, G. C. 2011, 

A&A, 531, A73 
Beccari, G., et al. 2010, ApJ, 720, 1108 
Bertin, E., & Amouts, S. 1996, A&A Supp., 117, 393 
Bessel, M.S., Castelli, E, & Plez, B. 1998, A&A, 333, 231 
Bopp, B.W., & Schmitz, M., 1978, PASP, 90, 531 
Boss, A. P. 1997, Science, 276, 1836 
Boss, A. P. 2003, ApJ, 599, 577 
Bouchy, E, et al. 2009, A&A, 496, 527 

Cai, K., Durisen, R. H, Michael, S., Boley, A. C, Mejfa, A. C, 
Pickett, M. K., & D'Alessio, P. 2006, ApJl, 636, L149 

Cai, K., Durisen, R. H, & Zhu, Z. 2009, Bulletin of the American 
Astronomical Society, 41, #439.01 

Calvet, N., Hartmann, L. & Strom, E. 2000, in "Protostars and 
Planets", eds V. Mannings, A. Boss, S. Russell (Tucson: Univer- 
sity of Arizona Press), 377 

Calvet, N., Muzerolle, J., Briceno, C, Hernandez, J., Hartmann, 
L., Saucedo, J. L., & Gordon, K. D. 2004, AJ, 128, 1294 

Clarke, C. & Pringle, J. 2006, MNRAS, 370, L10 

Cossins, P., Lodato, G., & Clarke, C. J. 2009, MNRAS, 393, 1 157 

Dahm, S. 2008, AJ, 136, 521 

Da Rio, N., Gouliermis, D. A., & Gennaro, M. 2010, ApJ, 723, 
166 

De Marchi, G., Panagia, N., & Romaniello, M. 2010, ApJ, 715, 1 
(Paper I) 

De Marchi, G., Panagia, N., Romaniello, M., Sabbi, E., Siri- 
anni, M., Prada Moroni, P. G., & DeglTnnocenti, S. 2011, 
larXiv: 1104.44941 (Paper H) 
De Marchi, G., et al. 2011, |arXiv:1106.2801| 
De Marchi, G., Panagia, N., & Sabbi, E. 2011, |arXiv:1106.5780l 
Dullemond, C. P., & Dominik, C. 2005, A&A, 434, 971 
Durisen, R. H., Boss, A. P., Mayer, L., Nelson, A. E, Quinn, T, & 

Rice, W. K. M. 2007, Protostars and Planets V, 607 
Elmegreen, B. G. 2000, ApJ, 530, 277 

Ercolano, B., Clarke, C. J., & Drake, J. J. 2009, ApJ, 699, 1639 
Ercolano, B., & Clarke, C.J., 2010, MNRAS, 402, 2735 
Ercolano, B., & Owen, J. E. 2010, MNRAS, 406, 1553 



Esteban, C, & Peimbert, M. 1995, Revista Mexicana de Astrono- 

mia y Astrofisica Conference Series, 3, 133 
Faesi, C, Covey, K. R., Gutermuth, R., et al. 201 1, Bulletin of the 

American Astronomical Society, 43, #340.32 
Fang, M., van Boekel, R., Wang, W., et al. 2009, A&A, 504, 461 
Fedele, D., van den Ancker, M.E., Henning, T, Jayawardhana, R., 

& Oliveira, J.M., 2010, A&A, 510, A72 
Fischer, D.A., & Valenti, J. 2005, ApJ, 622, 110 
Fitzpatrick, E. L., & Savage, B. D. 1984, ApJ, 279, 578 
Garmany, C. D. 1998, Origins, 148, 184 
Gonzalez, G. 1997, MNRAS, 285, 403 
Gorti, U„ & Hollenbach, D. 2009, ApJ, 690, 1539 
Gullbring, E., Hartmann, L., Briceno, C, & Calvet, N. 1998, ApJ, 

492, 323 

Haisch, K.E., Lada, E.A., & Lada, C.J. 2001, ApJ 
Hartigan, P., Kenyon, S.J., Hartmann, L., et al. 1991, ApJ, 382, 
617 

Hartmann, L., Calvet, P., Gullbring, E., D'Alessio, P. 1998, ApJ, 
495, 385 

Hartmann, L., 2003, ApJ, 585, 398 

Hartmann, L., D'Alessio, P., Calvet, N., & Muzerolle, J. 2006, 
ApJ, 648, 484 

Hartmann, L., Zhu, Z., & Calvet, N. 2011. larXiv:1106.3343l 
Harvey, P., Merin, B., Huard, T.L., et al. 2007, ApJ, 663, 1149 
Haywood, M. 2008, A&A, 482, 673 
Herbig, G. H. 1960, ApJs, 4, 337 
Herbst, W., 1986, PASP, 98, 1088 
Herczeg, G.J, & Hillenbrand, LA., 2008, ApJ 681, 594 
Hester, J.J., Scowen, PA., Sankrit, R., et al. 1996, AJ, 111, 2349 
Hollenbach, D.J., Yorke, H.W., & Johnstone, D. 2000, Protostars 

and Planets IV, 401 
Ida, S., & Lin, D. N. C. 2004, ApJ, 616, 567 
Isella, A., Carpenter, J.M., & Sargent, A.I., 2009, ApJ, 701, 260 
Jayawardhana, R., Coffey, J., Scholz, A., et al. 2006, ApJ 648, 

1206 

Jeffries, R. D. 2011, |arXiv:1102.4752| 
Johansen, A., & Youdin, A. 2007, ApJ, 662, 627 
Johansen, A., Youdin, A., & Mac Low, M.-M. 2009, ApJl, 704, 
L75 

Kenyon, S. J., & Hartmann, L. 1995, ApJS, 101, 117 
Kim, K. H, et al. 2009, ApJ, 700, 1017 
Kraus, A. L., & Hillenbrand, L. A. 2008, ApJl, 686, LI 1 1 
Linsky, J.L., Gagne, M., Mytyk, A., et al. 2007, ApJ, 654, 347 
Maeder, A., Grebel, E.K., & Mermilliod, J.-C. 1999, A&A, 346, 
459 

Marcy, G. W„ & Butler, R. P. 1998, Ann. Rev. A&A, 36, 57 

Marigo, P., Girardi, L., Bressan, A., et al. 2008, A&A, 482, 883 

Mayne, N. J., Naylor, T, Littlefair, S. P., Saunders, E. S., & Jef- 
fries, R. D. 2007, MNRAS, 375, 1220 

Mayor, M., & Queloz, D. 1995, Nature, 378, 355 

Mercer, E. P., Miller, J. M., Calvet, N., Hartmann, L., Hernandez, 
J., Sicilia-Aguilar, A., & Gutermuth, R. 2009, AJ, 138, 7 

Merin, B., Jorgensen, J., Spezzi, L., et al. 2008, ApJS, 177, 551 

Merin, B., et al. 2010, ApJ, 718, 1200 

Milone, A.P, Bedin, L.R., Piotto, G., & Anderson, J., 2009, A&A, 
497, 755 

Mohanty, S., Basri, G., & Jayawardhana, R. 2005, Astronomische 

Nachrichten, 326, 891 
Moor, A., Abraham, P., Juhasz, A., et al. 201 L larXiv: 1109.22991 
Muzerolle, J., Calvet, N., Briceno, C, Hartmann, L. & Hillen- 
brand, L. 2000, ApJ, 535, L47 



Mass accretion rates ofPMS stars in the LMC 21 



Muzerolle, J., Luhman, K., Briceo, C, Hartmann, L. & Calvet, N. 

2005, ApJ, 625, 906 
Natta, A., Testi, L., Muzerolle, J., Randich, S., et al. 2004, A&A, 

424, 603 

Natta, A., Testi, L., & Randich, S. 2006, A&A, 452, 245 
Nguyen, D.C., Scholz, A., van Kerkwijk, M.H., et al. 2009, ApJ, 
694, L153 

Owen, J. E., Ercolano, B., Clarke, C. J., & Alexander, R. D. 2010, 

MNRAS, 401, 1415 
Padgett, D. 1996, ApJ, 471, 847 

Panagia, N. 1991, in "New Views of the Magellanic Clouds", IAU 

Symp. 190, eds. Y.-H. Chu, N. Suntzeff, J. Hesser, D. Bohlender 

(San Francisco: ASP), 549 
Panagia, N., Gilmozzi, R.,Macchetto, E, Adorf, H.-M. & Kirsh- 

ner, R. P. 1991, ApJ, 380, L23 
Panagia, N., Romaniello, M., Scuderi, S. & Kirshner, R. 2000, 

ApJ, 539, 197 

Patel, N. A., Goldsmith, P. E, Heyer, M. H., Snell, R., & Pratap, 
P. 1998, IAU Colloq. 166: The Local Bubble and Beyond, 506, 
42 

Petrov, P. P., Gahm, G. E, Stempels, H. C, Walter, F. M., & Arte- 
menko, S. A. 201 l, |arXiv:1109.1266l 

Pinzon, G, de la Reza, R., Torres, C.A.O., Seifahrt, A., 
Neuhaiiser, R., Proceedings of the IAU Symposium 243 - Star- 
disk interaction in young stars, Grenoble, 21-25 May, 2007 

Rigliaco, E., Natta, A., Randich, S., Testi, L., & Biazzo, K. 2011, 
A&A, 525, A47 

Romaniello, M. 1998, PhD thesis, Scuola Normale Superiore, 
Pisa, Italy 

Romaniello, M., Panagia, N., Scuderi, S. & Kirshner, R. 2002, AJ, 

123, 9157 
Salpeter, E. E. 1955, ApJ, 121, 161 

Santos, N.C., Israelian, G, Mayor, M., Rebolo, R., & Udry, S. 

2003, A&A, 398, 363 
Santos, N. C, Mayor, M., Bouchy, E, Pepe, E, Queloz, D., & 

Udry, S. 2007, A&A, 474, 647 
Scuderi, S., Panagia, N., Gilmozzi, R., Challis, P. & Kirshner, R. 

1996, ApJ 465, 956 
Setiawan, J., et al. 2003, A&A, 398, L19 

Setiawan, J., Klement, R. J., Henning, T, Rix, H.-W., Rochau, B., 
Rodmann, J., & Schulze-Hartung, T. 2010, Science, 330, 1642 

Sicilia-Aguilar, A., Hartmann, L., Hernandez, J., Briceno, C, Cal- 
vet, N. 2005, AJ, 130, 188 

Sicilia-Aguilar, A., Hartmann, L.W., Furesz, G. et al. 2006, AJ, 
132, 2135 

Sicilia-Aguilar, A., Henning, T, & Hartmann, L.W., 2010, ApJ, 
710, 597 

Skrutskie, M.E, Cutri, R.M., Stiening, R. et al. 2006, AJ, 131, 
1163 

Spezzi, L., Alcala, J.M., Covino, E., et al. 2008, ApJ, 680, 1295 
Tognelli, E., Prada Moroni, P. G., & DeglTnnocenti, S. 2011, 
larXiv: 1107.23181 

Turner, N. J., Sano, T, & Dziourkevitch, N. 2007, ApJ, 659, 729 
Turner, N. J., & Sano, T. 2008, ApJl, 679, L131 
Vinko, J., Sarneczky, K, Balog, Z., et al., 2009, ApJ, 695, 619 
Yasui, C, Kobayashi, N., Tokunaga, A.T., Saito, M., & Tokoku, 
C, ApJ, 705, 54 

Waters, L. B. F. M., & Waelkens, C. 1998, Ann. Rev. A&A, 36, 
233 

Wegner, W. 2000, MNRAS, 319, 771 
White, R.J., & Basri, G. 2003, ApJ, 582, 1109 
White, R„ & Hillenbrand, L. 2004, ApJ, 616, 998 



Wilking, B. A., Gagne, M., & Allen, L. E. 2008, Handbook of Star 

Forming Regions, Volume II, 351 
Worden, S.P., & Peterson, B.M., 1976, ApJl, 206, L145 
Zaritsky, D., Harris, J., Thompson, I.B., & Grebel, E.K. 2004, AJ, 

128, 1606 

Zhu, Z., Hartmann, L., & Gammie, C. 2010, ApJ, 713, 1 143 

This paper has been typeset from a Tj5X/ ET]gX file prepared by the 
author. 



